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1.0  INTRODUCTION  AND  SUMMARY 


A  military  need  exists  for  a  microwave  solid  state  device  technology 
with  a  complementary  analog  signal  processing  capabilities  to  that 
developed  over  the  past  ten  years  in  surface  acoustic  waves  (SAW) 
technology.  In  this  technology,  non-recursive  transversal  filters  have 
been  realized  through  either  interdigital  transducers  and/or  periodic 
reflective  arrays  defined  on  the  piezoelectric  substrates  in  which  the 
surface  acoustic  waves  are  propagated.  Utilizing  this  technology  simple 
delay  lines,  complex  matched  filters,  high  Q  resonators  and  chirp  filters, 
have  been  realized  hi  i he  VHF,  UHF  frequency  range.  Unfortunately,  the 
extension  of  this  technology  to  the  microwave  region  ( f > 1  GHz)  has  been 
difficult.  At  these  frequencies  the  surface  acoustic  wave  length  is  less 
than  1pm  making  fabrication  difficult  and  the  propagation  loss  becomes 
excessive  (>100db/  usee  at  10  GHz.) 

This  final  report  for  the  third  year  of  contract  AFOSR  80-0264 
summarizes  the  results  of  a  fundamental  investigation  of  a  viable 
technological  realization  of  this  goal  through  magnetostatic  waves  MSW 
propagating  at  microwave  frequency  in  magnetically  biased,  liquid  phase 
epitaxial  films  of  yttrium  iron  garnet  (YIG)  grown  on  gadolinium  gallium 
garnet  (GGG).  This  technology  has  a  number  of  advantages;  low  loss  (>30db/ 
usee  at  xband),  tunable  by  bias  field  to  any  center  frequency  (l-20GHz), 
long  wavelengths  (10-lOOpm  at  any  frequency),  and  finally  a  well  understood 
and  characterized  wave  phenomena  (MSW)  is  utilized.  The  necessary 
non-recursive  transversal  filtering  is  achieved  through  multi-element 
transducer  and/or  the  interaction  of  the  MSW  with  periodic  structures  in  or 


This  report  summarizes  the  work  accomplished  in  the  last  year  toward 
achieving  MSW  transversal  filtering.  First  studies  of  MSSW  (Magnetostatic 
Surface  Wave)  transducers  composed  of  periodic  arrays  of  narrow  shorted 
microstrips  are  summarized.  These  studies  develop  an  empirical 
(experimentally  derived)  3  port  model  for  this  type  of  transducer.  Next, 
studies  of  a  variable  time  delay  device  based  on  a  cascaded  MSSW  forward 
wave  line  and  an  MSBVW  backward  wave  delay  line,  are  presented.  Finally, 
studies  of  MSSW  resonators  and  resonators  oscillators  are  summarized. 


2.0  THREE-PORT  SCATTERING  PARAMETER  MODELING  FOR  MSSW 


ARRAY  TRANSDUCERS 


This  research  presents  work  accomplished  to  model 
Magnetostatic  Surface  Wave  Array  Transducers.  This  work 
starts  by  deriving  a  3-port  scattering  parameter  model  by 
use  of  several  previous  models.  The  derivation  is 
presented  in  two  parts,  namely  the  wave  ports  and  the 
electric  ports.  The  results  are  shown  for  single  and 
multiple  bar  arrays  along  with  the  equivalent  experimental 
results.  The  data  given  by  the  derived  model  will  show 
that  this  model  is  only  good  for  up  to  two  bar  arrays. 
Next  an  empirical  (experimentally  derived)  3-port  model 
will  be  investigated.  The  layout  of  several  devices  which 
were  used  to  measure  the  3-port  scattering  parameters  for 


CHAPTER  II 


2.1  THEORETICAL  3-PORT  MODEL 


2.1.1  introduction 

This  section  outlines  the  steps  in  the  development  of 
a  theoretical  3-port  scattering  parameter  model  for  MSSW 
transducers.  A  computer  program  has  been  developed  based 
on  this  model.  Limitations  of  the  model  will  be  discussed 
as  well  as  recommendations  for  further  work. 


2.1.2  Derivation  and  Development 


Any  linear  n-port  electrical  network  may  be  analyzed 
in  several  ways.  Convenience  usually  dictates  wh-ich 
method  is  used.  For  example,  if  the  network  operates  open 
circuited,  impedance  parameters  are  best,  but  if  it  is 
operated  shorted  circuited,  admittance  parameters  are 
best.  Now,  when  the  network  operates  under  matched  load 
conditions,  scattering  parameters  are  usually  used. 
Scattering  parameters  form  a  transformation  between 
variables  which  are  linear  combinations  of  voltages  and 
currents  in  a  network.  They  are  very  convenient  for 
problems  which  involve  insertion  loss  and  matching 
networks.  [7] 


The  scattering  parameters  for  MSSW  have  been  derived 
using  techniques  similar  to  those  of  Hasnain  [6]  for  MSFVW 
or  Smith  [8]  for  SAW  Because  Wu's  f 5 ]  impedance  model 
based  on  the  Poynting  vector  is  used  in  deriving  the  S- 
parameters,  a  power  normalization  was  chosen.  The  3-port 
model  has  two  wave  ports  and  one  electric  port  as  shown  in 
Figure  2.1.  Since  MSSW  is  anisotropic  and  is 
characterized  by  two  surface  waves  traveling  on  the  two 
ferrite  surfaces  in  opposite  directions,  the  forward  wave 
will  be  defined  as,  "the  wave  traveling  from  port  1  to 
port  2".  A  Transmission  line  model  for  MSW  derived  by 
Snapka  [9]  will  be  used  to  find  wave  port  S-par ameter s. 
This  model  consists  of  forward  and  reverse  wave 
transmission  and  reflection  coefficients  which  can  be 
easily  derived  using  Kirchoff's  laws  for  voltage  and 
current.  The  electric  port  S-parameters  will  be  found  from 


Figure  2.2a:  Conditions  for  Finding  and  s 


Fiaure  2.2t:  Conditions  for  Finding  S  and  S 


where  Zp  and  ZR  are  the  forward  and  reverse  wave  impedance 
and  _1  F  and] R  are  forward  and  reverse  propagation 
constants  If  an  open  circuited  microstrip  is  assumed  as 
wave  loading,  we  can  calculate  S^,  S21  and  S22- 

match  terminating  port  2  in  its'  characteristic  impedance 
and  exciting  port  1  {see  Figure  2.2a),  we  can  calculate 
S 1 1  a  nd  S  2 \ ■ 


Where : 


S11 


(2.2) 


S21  * 


b2 

al 


V1  =  al\/z7  + 


(2.3) 


(2.4) 


(2.5) 


(2.6) 


And : 


S11 


A  -  ZFC  -  D  +  B/Zp 
—  Z  n  A 


B 


Zi?D  +  ZcZoC 


(2.7) 


( BC  -  AD) ( ZF  +  ZR) 
ZpZpC  —  ZpD  +  B  —  ZpA 


(2.8) 


Then  by  inverting  the  transmission  matrix,  matching 
terminating  port  1,  and  exciting  port  2  (see  Figure  2.2b.) 
S12  and  S22  can  be  found: 


(2.9) 


’22  = 

a2 


V2  =  32\Ar  +  b2\f*F 


(2.10) 


(2.11) 


(2.12) 


V1  *  t>lv/ZR 


(2.13) 


(2.14) 


And : 


~  ( +  ZR ) 


ZpZoC  “  ZpA  ~  ZpD  +  B 


(2.15) 


-'1 


*.  *.  \  .  *.  \ 


D  +  B/Zd  -  ZpC  -  A 


Z pZ pC  —  ZpA  -  ZpD  +  B 


(2.16) 


12  was  chosen  to  flow  in  the  opposite  direction  of  its1 
definition  for  cascading  purposes. 


2.1.4  Electr ic  Port  S-parameters 

The  electric  port  S-parameters  were  derived  by 
assuming  the  wave  ports  were  match  terminated  (see  Figure 
2.3)  and  that  the  electric  port  is  also  matched  by 
adjusting  the  thickness  of  the  microstrip.  Adjusting  the 
thickness  of  the  microstrip  allows  one  to  find  the  value 
of  conductor  resistance  that  will  give  intrinsically 
better  coupling,  hence  a  better  match  to  the  driving 
electromagnetic  field.  [4]  To  derive  the  electric  port  S- 
parameters,  Wu's  impedance  model  was  used  (see  Figure 
2.4.).  Wu  calculated  the  equivalent  circuit  impedances  by 
use  of  the  Poynting  vector.  In  his  model,  Rm^  Xm^,  Pm2, 
and  Xm2  are  the  equivalent  forward  wave  resistance, 
forward  wave  reactance,  reverse  wave  resistance,  and 
reverse  wave  reactance  respectively  of  the  MSSW.  Rc  and 
Xc  are  due  to  the  conduction  loss  in  the  microstrip,  with 
ports  1  and  2  matched,  where  R2  is  the  source  resistance, 
the  electric  port  parameters  are  found. 
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(2.34) 


At  present  there  is  no  direct  way  to  find  S31  or  S32  (wave 
port  to  electric  port)  so  we  assume  S23  *  S33  and  S32  * 
S33*  The  only  justification  for  this  assumption  is  that 
S23  and  S31  are  both  forward  waves  and  S32  and  S33  are 
both  reverse  waves.  However,  there  is  a  difference  in  the 
coupling  in  that  S23  and  S33  represent  electromagnetic  to 
magnetostatic  coupling  and  S32  and  S33  represent 
magnetostatic  co  electromagnetic  coupling.  Having  all  the 
S-parameters,  we  next  convert  them  to  Y-parameters.  The 
conversion  to  Y-parameters  is  done  by  assuming  the  S- 
parameters  are  terminated  in  the  characteristic  impedance 
of  the  unmetalized  region.  Since  the  conversion  from  S- 
to  Y-parameters  is  straight  forward,  it  will  not  be  shown 
here.  [6]  Finally,  cascading  of  the  Y-parameters  was  done 
in  order  to  model  arrays  and  achieve  the  overall  response 
of  the  coupling  from  the  input  transducer  array  to  the  YIG 
and  then  to  the  output  transducer  array  (see  Figure  2.5). 
These  eguations  were  also  straight  forward  since  all  they 
involved  was  writing  a  system  of  equations  and  solving  by 
Laplace's  expansion.  These  equations  from  the  S- 


parameters  on  down  to  the  response  function  were  derived 


and  programmed  in  Fortran  on  an  IBM  4341  and  are  shown  in 
Appendix  1.  The  dispersion  relationship  for  MSSW  utilized 
is  that  of  Brinlee.  (1/3)] 

2.1.5  Results  and  Problems  Encountered 

The  simplest  case  of  a  single  bar  50  micrometers  wide 
transducer  was  first  tried.  The  computer  program  appeared 
to  yield  the  correct  frequency  response  as  that  found  from 
experiment  (see  Figures  2.6a  and  2.6b).  Then  cascaded 
microstrips  both  2-  and  4-bar  arrays  were  tried  and  the 
program  did  not  yield  the  correct  frequency  response. 
Figures  2.7  and  2.8  show  the  experimental  and  theoretical 
frequency  response  of  a  2-bar  and  4-bar  transducer  array 
pairs.  The  frequency  responses  of  these  devices  had  gain 
at  some  frequencies  which  is  not  physically  possible  for 
passive  devices.  Conservation  of  energy  was  then 
considered  and  the  squared  magnitudes  were  output  by  the 
program.  Conservation  of  energy  for  a  lossless  3-port 
says : 


S11 

2  + 

S31 

2  + 

S  2 1 

S, 

2  + 

t 

s  1  2 

2  + 

S32 

S  3  3 

2  + 

S1 3 

2  + 

S2  3 

(2.35) 

(2.36) 

(2.37) 
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Magnitude  in  PB 


Since  this  is  not  the  lossless  case,  the  sum  of  the 
squared  magnitudes  should  be  less  than  one.  As  things 
turned  out,  all  the  squared  magnitudes  were  indeed  less 
than  one  but  the  sums  were  not.  The  first  two  equations 
2.35  and  2.36  were  violated  while  equation  2-37  was  not 
because  it  came  directly  from  Wu's  impedance  model  derived 
from  the  Poynting  vector.  The  equations  2-35  and  2.36 
were  violated  because  they  came  from  the  open  circuited 
microstrip  transmission  line  model.  The  transmission  line 
model  behaves  very  much  like  an  ideal  transmission  line  in 
that  *  S22  *  0  an<3  S21  *  s12-  From  a  physical  point 
of  view  this  is  not  possible  since  there  must  be  some 
coupling  to  the  electric  port.  Some  other  approaches  were 
tried*  in  order  to  overcome  these  problems  including 
forcing  equations  2.35  and  2.36  to  obey  conservation  of 
energy  by  adjusting  the  magnitude  of  S^2  and  S21  while 
leaving  the  phase  alone.  This  was  done  by  allowing 
equation  2.37  to  set  the  total  energy  and  then  solving  for 
the  new  magnitudes.  After  all  this  was  J..!  -,  gain  still 
existed  at  certain  frequencies.  The  Y-parameters  were 
then  converted  to  Z-parameters  and  some  Z-parameters  had 
negative  real  parts.  The  problem  seems  to  be  in  the  phase 
of  the  S-parameters  or  Y-parameters.  The  possibility  of 
the  computer  not  handling  the  complex  numbers  riqht 
exists,  however  the  real  problem  is  believed  to  be  in 
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finding  accurate  values  for  S  3  j  and  S32-  The  final 
solution  will  be  to  make  Poyntina  vector  calculations  for 
a  loaded  microstrip  and  "that"  solution  is  still  pendinq 
at  this  time.  Since  solving  the  loaded  microstrip  problem 
is  very  difficult,  in  the  following  chapters  an  empirical 
approach  will  be  tried.  A  device  structure  will  be 
developed  and  made  in  order  to  measure  the  S-parameters  of 
a  transducer  array. 


CHAPTER  III 


2.2  DEVICE  DESIGN  AND  MANUFACTURE 

.2.1  Introduction 

The  simplest  and  easist  built  transducer  for 
conversion  of  electromagnetic  to  magnetostatic  energy  is 
one  that  has  the  microstrip  fabricated  on  a  dielectric 
substrate  with  the  ferrite  flipped  over  on  it.  All 
transducer  designs  will  be  based  on  this  flipped 
configuration  as  shown  in  Figure  2.9.  Efficient 
conversion  of  enecgy  is  obtained  when  the  length  of  the 
microstrip  is  small  compared  to  the  electromagnetic 
wavelength,  and  when  the  width  is  small  compared  to  the 
magnetostatic  wavelength.  If  a  microstrip  is  excited  by  a 
current  and  the  ferrite  is  biased  by  a  magnetic  field,  two 
magnetostatic  waves  are  launched  by  the  microstrip  in 
opposite  directions  perpendicular  to  the  axis  of  the 
microstrip  (Figure  2.10). 

.2.2  Device  Design 

The  transducer  design  for  measuring  the  3-port  S- 
parameters  of  a  microstrip  transducer  array  is  very 
straight  forward.  Two  measuring  single  microstrip 
transducers  are  placed  symmetrically  on  both  sides  of  the 


microstrip  transducer  array  to  be  measured.  The  sinqle 
microstrip  transducers  are  placed  at  a  distance  greater 
than  that  of  a  magnetostatic  wavelength,  and  at  the  same 
time,  not  equal  to  any  integer  multiple  of  a  quarter  of 
that  same  wavelength.  All  total  sixteen  different  devices 
were  made  including  single,  two,  four,  and  eight  bar 
microstrip  arrays  with  widths  of  five,  ten,  twenty,  and 
fifty  micrometers.  The  measuring  microstrip  transducer 
width  was  chosen  as  five  micrometers  because  it  was  the 
narrowest  line  that  could  be  fabricated  consistently. 
Narrower  transducers  should  have  a  tendency  to  be  more 
decoupled  from  the  microstrip  transducer  array  under 
measut ement . 

In  the  single  and  two  bar  microstrip  arrays,  the 
measuring  transducers  were  2200  micrometers  or  seven  and 
one-third  magnetostatic  wavelengths  away  from  the  array; 
while  in  the  four  and  eight  bar  arrays,  the  measuring 
transducers  were  1600  micrometers  a  way  or  five  and  one- 
third  wavelengths.  These  distances  were  chosen  to  meet 
the  constraints  of  the  above  and  for  ease  of  fabrication. 
The  center-to-center  spacing  of  the  input  and  output 
arrays  was  chosen  to  be  one  centimeter  for  convenience, 
and  the  microstrip  length  was  three  millimeters  or  one- 
tenth  of  an  electromagnetic  wavelength.  The  transmission 
lines  connecting  all  transducer  arrays  to  the  outside 


world  were  50  Ohms  for  measurement  purposes.  Fiaure  2.11 
shows  the  layout  of  a  typical  device.  All  of  the  devices 
were  then  fabricated  and  tested. 

2.2.3  Dev  ice  Fabrication 

After  the  transducer  desian  was  completed, 
fabrication  of  the  devices  was  next.  Special  care  had  to 
be  taken  in  making  these  devices  since  lengthy  five 
micrometer  wide  lines  are  not  easy  to  achieve  with  a 
single  reduction  photographic  Drocess.  A  Haag-Streit  Bern 
coord  inatograph  was  used  to  cut  masks  for  each  device. 
This  coordinatograph  could  cut  all  lines  except  the  five 
micrometer  wide  lines  because  its'  resolution  is  about 
ten  micrometers.  To  replace  the  five  micrometer  lines, 
fine  copper  wire  of  the  appropriate  width  was  used.  This 
wire  was  straightened  and  fixed  to  the  mask  by  tape  at  the 
ends  and  clear  silicon  grease  in  the  middle.  Each  mask 
was  shot  with  a  25.4  to  1  reduction  camera  where  focusing 
was  tried  prior  to  the  reduction  process,  but  was 
abandoned  since  it  could  only  be  achieved  over  a  small 
region  of  the  mask.  The  lack  of  a  good  focus  caused  many 
problems-  which  were  overcome  by  under-developing  the 
photographic  reduction  by  about  fifty  percent.  The 
reductions  of  each  mask  were  carefully  examined  for 
problems.  Each  line  was  checked  for  unifqrmity,  width, 
and  spacing  with  a  Lietz  microscope.  An  actual 


reproduction  of  a  reduced  mask  can  be  seen  in  Figure  2.12 
After  the  photo-reduction  process  came  the  photo-resist 
process . 

KTI  positive  photo-resist  with  a  viscosity  of  90 
centistokes  @25°  C  was  spun  onto  a  substrate  at  5000  rpm 
for  40  sec.  yielding  a  thickness  of  about  3  micrometers. 
The  photo-resist  was  then  prebaked  for  about  15  minutes  at 
70°  C.  After  the  prebake,  the  substrates  were  exposed 
with  a  Colbilt  CA-400  mask  aligner.  The  exposure  time 
depended  on  the  quality  of  each  mask.  The  substrates  were 
then  post  baked  overnight  for  a  hard  bake  to  prevent 
widing  of  the  lines  in  the  up-plating  process.  The  10 
mil.  alumina  substrates  have  about  100  Angstroms  of  chrome 
and  about  700  Angstroms  of  gold  evaporated  on  both  sides. 
Before  the  evaporation  process,  the  substrates  were 
cleaned  with  a  one  to  one  solution  of  sulfuric  acid  and 
hydrogen  peroxide  and  baked  dry.  After  the  substrates 
were  cleaned  and  dried,  they  were  placed  in  the 
evaporation  system  and  futher  dried  with  the  substrate 
heater 

The  up-plating  process  consisted  of  Englehard  ECF-63 
noncyanide  gold  plating  solution,  a  hot  plate  stirrer,  and 
an  HP  model  721A  power  supply.  The  substrates  (the 
cathode)  were  clamped  in  parallel  each  with  a  alligator 
clip  in  series  with  a  IK  Ohm  resistor,  a- circular  platinum 


electrode  was  submerged  in  the  plating  solution  and  served 
as  the  anode.  The  solution  required  that  no  glass  be  used 
so  a  nalgene  beaker  served  as  the  container  for  the 
solution  and  the  electrodes  (see  Figure  2. 13}  The 
solution  was  first  filtered  thru  a  Gelman  syringe  filter 
and  then  heated  to  45-50  °  C.  A  current  of  10  m a  per 
substrate  for  30  min.  yielded  about  3  micrometers  of  gold. 
A  problem  was  encountered  with  the  up-plating,  in  that, 
the  gold  was  not  sticking  very  well.  This  problem  was 
thought  to  be  caused  by  unclean  substrates  and  was 
resolved  by  reversing  the  current  for  about  1  minute 
before  beginning  the  up-plating  process  and  every  10 
minutes  thereafter.  When  the  up-plating  was  complete,  the 
solution  was  refiltered  back  into  the  original  container 
and  the  hard  baked  photo-resist  was  stripped  off  with  a 
one-to-one  solution  of  sulfuric  acid  and  hydrogen 
peroxide.  The  substrates  were  etched  in  a  premixed 
solution  of  125  grams  of  KI,  40  grams  of  I2 •  and  500 
milliliters  of  water  to  remove  the  excess  gold  and  then 
etched  in  another  premixed  solution  of  50  qrams  of  NaOH, 
40  grams  of  {^FefCNJg,  and  180  milliliters  of  water  to 
remove  the  chrome.  After  careful  examination  under  a 
microscope  for  correct  width  and  line  continuity,  the 
substrates  were  mounted  on  an  aluminum  block  and  tested. 


The  3-port  S-parameters  represented  by  eauation  3.1 


are  to  be  measured. 


*■* 

— 

— 

-  — 

bl 

311 

3 1 2 

S13 

al 

b2 

= 

S21 

s  22 

S23 

a2 

(3.1) 

b3 

S31 

S32 

S3  3 

a3 

Previous  mathmatical  models  [6]  assumed  that  the  two  wave 
ports  (Fig.  2.14  )  were  match  terminated  making  the 
reflections  coefficents  zero.  Since  all  indications  from 
the  measurements  show  that  this  condition  holds,  and 
because  the  wave  port  reflection  coefficients  are  very 
hard  to  measure,  S^  and  S22  be  assumed  equal  to 
zero.  The  validity  of  this  assumption  as  well  as  other 
assumptions  will  be  examined  later.  S-parameters  are 
similar  to  Y-  and  Z-parameters  in  that  in  order  to  measure 
these  parameters,  one  must  set  something  equal  to  zero. 
For  example,  to  measure  833,  a2  and  83  must  be  zero. 
Consequentially,  the  ratio  of  b3  to  a^  will  be  S33-  Since 
excitation  of  a  microstrip  in  the  presence  of  a 
magnetically  biased  ferrite  launches  waves  in  two 
directions,  one  must  develop  techniques  for  eliminating 
the  other  wave  in  order  to  measure  the  S-parameters.  Due 
to  the  fact  that  the  reverse  wave  is  more  difficult  to 
couple  to,  it  generally  contains  much  less  power  than  the 


analy2er,  which  has  the  ability  to  measure  reflection  and 
transmission  coefficients  directly. 

2.2.5  Conclusion 

This  chapter  showed  how  each  device  was  made  and 
measured.  All  the  measurements  taken  were  stored  on 
cassett  tape  and  later  transferred  to  an  IBM  4341  computer 
and  stored  on  disk.  The  following  chapter  studies  this 
data  and  examines  an  empirical  model  derived  from  it. 


CHAPTER  IV 


2.3  EMPIRICAL  3-PORT  MODEL 

2.3.1  Introduction 

Measuring  S-parameters  of  a  transducer  with  a  similar 
transducer  at  microwave  frequencies  introduces  errors. 
These  errors  can  be  corrected  by  introducing  correction 
factors.  The  errors  result  from  transducer  conversion 
loss,  the  propagation  path  between  the  measurina 
transducer  and  the  transducer  array  being  measured,  and 
the  transmission  line  between  the  network  analyzer 
reference  plane  and  the  transducer  reference  plane.  These 
correction  factors  will  be  derived  and  the  resulting 
measurements  will  be  corrected  and  placed  in  the  S- 
parameter  model  derived  in  Chapter  2. 

2.3.2  Correction  Factors 

The  transducer  conversion  loss  is  assumed  to  be  a 
real  constant  loss  and  is  accounted  for  by  adding  the 
assumed  values  of  5  dB  per  transducer  for  forward  wave 
coupling  and  10  dB  for  reverse  wave  coupling.  The  phase 
difference  between  the  reference  planes  was  measured  with 
the  network  analyzer  by  taking  an  equivalent  section  of 
transmission  line  and  shorting  the  end.  The  measured 


forward  wave,  therefore  a2  can  be  neglected  in  the 
measurement  of  S2j,  S23,  333,  and  S33-  However,  in 
measuring  the  reverse  wave  S- pa r am e t e r s ,  S1?,  S32,  and 
S13,  a^  cannot  be  neglected.  This  problem  can  be 
eliminated  by  cleaving  the  ferrite  crystal  at  one  end  and 
placing  the  cleaved  edge  directly  on  the  transmitting 
transducer  as  shown  in  Figure  2.15  This  forces  all  the 
energy  in  one  direction  eliminating  any  wave  that  might 
propagate  in  the  other  direction.  To  kill  secondary 
reflected  waves  created  by  unterminated  transducers,  50 
Ohm  loads  were  connected  to  each  transducer  not  currently 
being  used. 

Another  problem  associated  with  measuring  the  S- 
parameters  is  ripple  in  the  passband  due  to  reflections 
from  the  ends  of  the  ferrite  crystal.  This  problem  was 
partially  solved  by  terminating  the  end  of  the  crystal  not 
cleaved.  The  termination  was  achieved  by  lapping  the  end 
at  about  a  1°  angle  as  shown  in  Figure  2.15  Terminating 
the  ends  by  lapping  works  exceptionally  well  because  the 
thinner  the  ferrite  crystal  the  longer  the  delay 
associated  with  the  magnetostatic  wave.  Therefore,  any 
wave  traveling  to  the  end  is  attenuated  before  it  has  time 
to  be  reflected  back.  The  results  of  lapping  are  very 
good  as  can  been  seen  from  Figure  2.16  All  measurements 
were  taken  with  a  Hewlett-Packard  8<*09  automatic  network 


Figure  2.16  Unlapped  Crystal  Response 


angle  was  subtracted  from  130°  and  the  result  divided  by 
two,  yielding  the  correct  electrical  length  of  the 


transmission  line.  The  propagation  path  error  was 
corrected  by  solving  the  dispersion  relation  for  the 
propagation  constant  and  assuming  exponential  propagation- 
The  propagation  correction  factors  were  derived  as 
follows:  (refer  to  Fig.  2.14) 
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where  the  subscripts  "  m",  "  c " ,  " f " ,  and  "r"  represent  the 
measured  and  corrected  S-parameters,  and  the  forward  and 
reverse  waves  respectively. 

The  final  S-parameters  with  all  the  corrections  made 
are  as  follows: 


Sn*0  (4.13) 

S22  =  0  (4.14) 


s12=  (  s12m  +20  dB)*exp(2Crd)*exp[ j(/S12m+ 

(180-e1)/2  +  (180-©2)/2)]  (4.15) 

s21=  (  s21m  +  dB )  *exp(  2^1f  d )  *exp  [  j  ( /S2lm+ 


(180-©1)/2  +  (180-92)/2)]  (4.16) 

s13=  (  s13m  +1°  dB)*exp(S‘rd)*exp[  j(/S13m+ 

(180-91)/2  +  ( 180-©3 ) / 2  )  )  (4.17) 

s31=  (  S31m  +5  dB)*exp(I  fd)*exp[ j(/S31ro+ 

(  180-91)/2  +  (  1 8  0-©3 ) /2 ) ]  (4.18) 

^23=  (  s23m  +  5  dB ) *exp( Jfd ) *exp [ j ( /S23m+ 

( 1 80-©2 ) /2  +  ( 1 8  0-©3 ) /2 ) ]  (4.19) 

S32=  f  s 3 2m  +1°  dB)  *exp(/'rd  )*exp[  j  (/S32m+ 

(180-©2)/2  +  ( 1 80 -©3 ) /2 ) ]  (4.20) 


s33=  s33  exp[j(/S33m  +  ( 180-©3 ) ] 


(4.21) 


These  equations  were  programmed  into  an  IBM  4341  computer 
where  the  data  had  been  previously  stored. 

2.3.3  The  Measurement  System 

The  Hewlett-Packard  8409  Automatic  Network  Analyzer 
provides  reliable  accurate  measurements  of  a  2-port 
network.  Control  of  the  network  analyzer  was  done  by 
using  an  HP-9845B  Desk  Top  Computer  together  with  an 
Accuracy  Enhancement  Program.  The  Accuracy  Enhancement 
Program  was  written  by  Hewlett-Packard  and  has  the  ability 
to  automatically  measure  the  reflection  and  transmission 
coefficients  over  any  frequency  range  from  100  MHz.  to  18 
GHz.  and  then  take  the  measured  data  and  correct  it  using 
either  an  8-term  or  12-term  error  model.  All  measurements 
taken  for  this  thesis  used  the  12-term  error  model. 
Before  any  measurements  were  made,  the  system  had  to  be 
calibrated.  The  Accuracy  Enhancement  Program  provides  a 
calibration  routine  which  sets  the  reference  magnitude 
levels  and  the  phase  reference  plane  using  shorts,  opens, 
and  loads.  Since  the  device  being  measured  is  a  3-port 
device,  a  multi-port  measurement  system  would  be  more 
appropriate.  Such  systems  have  been  developed  but  not 
available  for  these  measurements.  The  limitations  of  the 
2-port  system  was  not  easily  overcome.  Each  time  a 
measurement  was  taken,  special  care  had  to  be  used  in 
lining  up  the  magnetically  tunable  passband.  For  all 
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measurements,  the  frequency  span  was  from  2.4  to  3.6  GHz. 
with  the  passband  centered  at  3  GHz.  The  data  was  stored 
on  tape  after  each  measurement  and  later  transferred  to 
the  IBM  for  processing.  The  device  under  test  was  placed 
on  a  galvanometer  for  proper  alignment  in  a*Varian  9-inch 
magnet . 

2.3.4  Measured  and  Corrected  Results 

The  first  measurements  were  on  a  single  microstrip, 
50  micrometer  wide  transducer.  The  results  are  shown  in 
Figure  2.18  (Note  the  similarity  in  S32,  S23  and  S13, 
S 3 2 - )  When  these  results  were  transferred  to  the  IBM,  S32 
was  set  equal  to  S23;  and  S  ^  3  was  set  equal  to  S32  to 
simplify  the  alignment  problems. 

The  following  pages  are  the  measured  results  of 
several  selected  devices.  Only  a  few  devices  were  chosen 
due  to  the  large  amount  of  data.  This  data  is  included  in 
Figures  2.18  thru  2.52  As  one  will  note,  there  is  a  small 
problem  with  alianment  of  the  measurements.  This  problem 
is  due  to  reversing  the  magnetic  field  and  is  somewhat 
overcome  by  shifting  the  frequency  of  the  measured  data 
before  using  it  in  the  model.  Figures  2.53  thru  2.87  shows 
the  corrected  data  of  the  devices  shown  in  Figures  2.18 
thru  2.52  As  one  can  see  from  comparing  the  corrected 
results  to  the  measured  results,  only  the  .angle  has 
actually  changed.  The  only  change  in  magnitude  that  can 
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be  detected  is  that  it  is  shifted  up  has  not  changed  but 
has  only  shifted  up. 

1.5  Empirical  Model  Results 

Figures  2.18  thru  2.52  show  the  frequency  response 
calculated  by  substituting  the  corrected  S-par ameter s  for 
the  theoretical  S-parameters  in  the  computer  program  of 
Appendix  1.  This  leads  to  an  empirical  model.  When  these 
results  are  compared  with  previous  experimental  data 
(Figures  2.93  thru  2.97  there  is  not  much  difference  in 
the  resonance  notch  structures  or  the  insertion  loss. 
This  indicates  that  the  model  is  good  at  least  to  the 
first  order.  To  exploit  the  model  more  fully,  one  must 
try  to  calculate  multiple  bar  responses.  Figures  2.98  and 
2.99  show  the  results  of  using  the  single  bar  50um  data  to 
predict  the  frequency  response  of  two  and  four  bar  arrays. 
The  calculated  two  bar  frequency  response  when  compared  to 
the  experimental  two  bar  frequency  response  of  Figure  2.100 
looks  good;  however,  the  calculated  four  bar  frequency 
response  does  not  match  the  experimental  four  bar 
frequency  response  of  Figure  2.10-1.  The  difference  is  in 
the  resonance  notch  structure  and  in  the  insertion  loss. 
The  resonance  structure  of  the  2-bar  array  looks  very 
close.  The  only  problem  is  in  the  roll  off  rates  of  each 
lobe.  This  indicates  that  the  phase  change  associated 
with  the  resonances  is  not  sharp  enough.  The  problem  may 
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be  helped  by  more  accurate  phase  combining  of  the  data. 
There  seems  to  be  many  small  problems  associated  with  the 
empirical  model.  The  most  pronounced  problem  is  probably 
measurement  errors  arising  from  misalignment  of  the 
different  measurements  on  each  device  and  the  inherent 
accuracy  of  the  network  analyzer.  Another  error  occurs 
from  the  fact  that  this  is  primarily  a  superoosition  model 
and  that  there  is  no  accounting  for  mutual  coupling 
between  the  microstrips.  These  errors  as  well  as  a  few 
other  problems  will  be  examined  in  the  concluding  chapter. 


CHAPTEF  V 


2.4  CONCLUSION 


2.4.1  introduction 

In  this  chapter,  a  final  look  will  be  taken  at  each 
model  and  their  results.  The  good  points  as  well  as  the 
weak  points  will  be  examined  for  both  models.  Finally, 
recommendations  for  futher  study  will  be  made. 

2.4.2  The  Theoretical  Model 

The  theoretical  model  of  Chapter  2  is  a  good  model 
for  simple  MSSW  structur.es.  Problems  arise,  however,  when 
the  structures  become  more  complex.  This  bre^k  down  is 
simply  due  to  the  assumptions  of  S23  =  S3^,  si3  =  s32'  an<^ 
that  the  transmission  line  model  can  be  used  to  find  the 
wave  port  parameters  directly.  These  assumptions  are 
necessary  due  to  the  very  complex  nature  of  MSSW  and  its' 
nonreciprocal  properties.  The  model  approach  of  using 
scattering  parameters  is  very  good  and  is  probably  the 
easiest  way  to  go  as  can  be  seen  from  the  empirical  model. 

2.4.3  The  Empir  ical  Model 

The  empirical  model  of  Chapter  4  is  an  excellent 
model  for  MSSW.  This  model  uses  the  same  principles  as 


the  theoretical  model.  The  only  difference  is  that  the 
empirical  model  uses  measured  S-par ameters ,  where  the 
theoretical  model  uses  mathematically  derived  parameters. 
The  empirical  model  is  not  without  its'  problems.  and 
S22  are  assumed  to  be  zero  because  they  could  not  be 
measured.  This  may  be  nearly  true  for  single  microstrip 
transducers,  but,  since  the  transducer  arrays  have 
interference  and  resonators  can  be  construced,  the 
condition  of  sn*s22=^  cannot  be  absolutely  true.  Most  of 
the  other  measured  S-parameters  are  as  expected  except  for 
the  appearance  that  S13*S32  and  S31=S23-  All  the  devices 
measured  (see  Figures  4.1  thru  4.5)  display  these 
characteristics.  The  reason  the  meaurements  show  this  is 
because  of  inadequate  error  correcting.  The  coupling  loss 
added  by  the  measuring  transducers  can  be  accounted  for, 
but  there  is  no  possible  way  to  separate  the  phase  change 
that  takes  place  in  going  from  electromagnetic  waves  to 
magnetostatic  waves  and  from  magnetostatic  waves  to 
electromagnetic  waves.  The  five  micrometer  transducers 
hopefully  are  decoupled,  but  there  will  still  be  some 
phase  change-  Another  error  in  this  model  is  that  it  does 
not  account  for  mutual  interaction  and  coupling  between 
the  microstrips  in  each  array.  Although  this  model  does 
have  some  problems,  it  is  still  able  to  yield  the  correct 
response  (see  Figures  2.88  thru  2.92  '  -  Difficulties  do 


arise  in  trying  to  cascade  the  single  microstrip  50 
micrometer  results  as  can  be  seen  in  Fiqures  2.98  and 
2.99.  The  two  bar  results  are  not  bad,  but  the  four  bar 
results  are.  The  reason  is  that  there  are  many  errors 
including  alignment  of  the  different  measurements  on  each 
device,  and  in  the  accuracy  of  the  automatic  network 
analyzer.  The  automatic  network  analyzer  is  a  good 
accurate  system  but  it  is  not  perfect.  The  errors  that  do 
occur  are  propagated  and  greatly  magnified  when  the 
microstrips  are  cascaded. 

2.4.4  Recommendations  for  Further  Study 

For  the  empirical  model,  there  is  probably  no  way  to 
improve  on  the  measured  data  unless  a  more  accurate 
measurement  system  is  found-  However,  study  of  the 
measurement  techniques  and  development  of  new  techniques 
might  possibly  yield  better  results.  The  real  solution  is 
to  address  the  theoretical  model  and  the  Droblems 
associated  with  it. 
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APPENDIX  I 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


COMPUTER  PROGRAM  FOR  CALCULATING  THE  THREE  PORT 
S-PARAMETERS  FOR  MSSW  TRANDUCERS.  THIS  PROGRAMS 
USE  THE  DISPERSION  PROGRAM  OF  BRINLEE  AND  THE 
RADIATION  IMPEDANCE  CALCULATED  FROM  WU ' S  MODEL. 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


DIMENSION  FRQ( 152 ) , WAV ( 152 ) ,PLOSS( 152 ) ,RGF1 ( 152 ) , 
$AGFl (152) 

DIMENSION  El ( 5 ) ,E2(5) ,D1(5) ,D2(5) ,SC1(5) ,SC2(5) 
DIMENSION  SI (25)  ,S2(25) ,AGF2 ( 152 ) ,AGR2 ( 152 ) 
DIMENSION  RGF2(152) , RGR2 ( 1 52 ) , RGR1 ( 1 52 ) ,AGR1(152) 
DIMENSION  RMF( 152 ) ,XMF<152) ,RMR(152) ,XMR(152) 
COMPLEX* 8  ARR1( 3,3) , ARR2 ( 3 , 3 ) , Al ( 3 , 3 ) ,A2( 3,3) , 

$Rl ( 3 , 3 ) , R2 ( 3 , 3  ) 

1, PP( 3,3) ,YSB( 3,3) , Z , VFl , VR1 , TZ1 , PPP ( 3,3) 

COMPLEX* 8  ZIN,CSINH,CCOSH, CTANH , JJ , TZ2 
COMPLEX*8  ZF2,ZR2,ZF1,SSB(3,3) ,ZRl 

2 ,  AM , BM , CM, DM, GAMMA, KFl , KRl ,KF2,KR2,W0 
COMPLEX* 8  GF1 ( 1 52 ) ,GR1(152) ,GF2(152) ,GR2(152)  , 

$Ul ,U2 , GAMMAF , GAMMAR 

C  INITIALIZATION  OF  VARIBLES 

DATA  P,KT,KAT, LL,M/1 .0E-05, 1 ,0 ,0 ,59/ 

DATA  El(l) , El ( 2 ) ,E1(3) , El v  4 ) ,E1(5)/10. 0,1. 0,1.0, 
$16.0,16.0/ 

DATA  E2(l) , E 2 ( 2 )  ,E2(3) , E2 ( 4 ) ,E2(5)/0. 0,0. 0,0.0, 
$0.0, 0.0/ 

DATA  SC 2(1)  ,  SC2 ( 2  )  , SC 2(3) ,SC2(4) , SC  2(5) /0. 0,0.0, 
$0.0, 0.0, 0.0/ 

DATA  IF, IFF, MIFF, I HEAD/0 ,1,0,0/ 

CS INH (Z)*(CEXP(Z)-CEXP(-Z) )/2. 

CCOSH( Z ) = ( CEXP ( Z)  CEXP(-Z) )/2. 

CTANH ( Z ) -CS INH ( Z ) /CCOSH ( Z ) 

DATA  EPS,NSIG,KN, NGS , N , ITMAX/1 . E-5 , 5 , 0 , 1  ,  1 , 50  / 
DATA  ZC,RS,DIFP/90  -  ,  50  .  ,0  ./ 


PI-ACOS ( -1  - ) 

U0-PI*4.E-7 
GYRO* 2. 8 
JJ»(0.0,1.0) 

READ (5,*)  D0,TU,T0,H0,SM,DHM,NSAMP 
READ( 5 ,  * )  T,B0,RLGH  OIST,Nl,N2 
READ (5,*)(S1(I),I*.'  Nl) 

READ(5,*)(S2(I) ,1-1, N2) 

READ  (5,*) (Dl( I) ,1=1,5) 

READ(5,*)  ( SC 1(1) ,1*1,5) 

DO  101  1*1,5 
101  D2 ( I ) =D1 ( I ) 

T1=B0 

DO  18  MM=1 , Nl 
18  Tl-Tl+Sl ( MM ) 

T2-B0 

DO  28  MM*1,N2 
28  T2-T2+S2 ( MM ) 

A0*DIST*1 .E3-(Tl+T2)/2 . 

WRITE( 6,66)  DO ,TU ,T0 ,H0 ,SM,DHM 
WRITE (6,67)  T , BO , RLGH , AO 
WRITE (6,68)  N1,(S1(I) ,1-1, Nl) 

WRITE (6,69)  N2,(S2(I) ,I*1,N2) 

KF2-CMPLX ( 0 . 393423885E-05 , 0  - 125572526E-03 ) 

KR2=CMPLX(0 .961030200E-05,0.307895565E-03) 

C  COMPUTE  PASSBAND 

FPM-1760 - 
HI-375. 

GA-2.8 

GM=GA*FPM 

OH-HI/FPM 

FRQ3-SQRT (OH* ( OH+1 . ) ) *GM 
FRQ4* ( OH+O -  5 ) *GM 
FRQ1=FRQ3*1 . 01 
FRQ2=FRQ4*0 .99 
QF-0 . 5* ( FRQ2+FRQ1 ) 

DF*(QF-FRQl)/100 -0 

C  FREQUENCY  SAMPLING  152  POINTS  FOR  DISPERSION  RELATION 

FRQ( 1 ) -FRQ3+0 . 75* ( FRQ 1- FRQ 3 ) 

FRQ ( 2 ) -FRQl+0 . 5*DF 
DO  102  1=1,49 
RI=  I 
J  =  I+2 

102  FRQ (J)=FRQ(2)+RI*DF*2-0 
FRQ( 52 ) =QF+DF*0 . 25 
DO  103  1=1,99 


3  FRQ( J)*FRQ( 52 )+RI*DF 

FRQ( 152 ) *FRQ2+0 -25*(FRQ4-FRQ2) 

FORMAT ( ’  FRQ'  ,  E16.8) 

WM-GYRO*SM 
GH»GYRO*HO 
DHM2*DHM/2 .0 
STARTF=FRQ{  D/1000 .0 
STOP F * FRQ ( 152)/1000 .0 
UM=PI*4 -0E-09 

BEGIN  FREQUENCY  SCAN 
NSAMP=NSAMP-1 
DO  10  M=1,NSAMP 
F=*FRQ(  M) 

PRSNTF=F/1000 .0 

DH=DHM2* (3000-0/F+F/ 3000.0  ) 

WO»CMPLX(GH, (GYRO*DH/2 .0 ) ) 

Ul=*l  ,0+W0*WM/(W0**2-F**2  ) 

U2=WM*F/( W0**2-F**2) 

KFl»JJ*CLOG( ( (Ul-U2-1-0)*(U1+U2) )/( ( U1+U2+1 . 0 ) * ( 

$Ul-U2) ) )/2.0/Dl(4) 

KF1=CMPLX( ( 3  - 0*REAL ( KFl ) ) , AIMAG ( KF1 ) ) 

KRl*JJ*CLOG( ( (U1+O2-1-0) *(01-02) )/( (U1-U2+1 .0 > * ( 

SU1+U2) ) )/2.0/Dl(4) 

KRl*CMPLX ( ( 2 . 0*REAL ( KRl ) ) r AIMAG ( KR1 ) ) 

CALL  DSPRSN( F,U1,U2,E1,D1,SC1, EPS , GAMMAR, GAMMAF , 

$KF 1 , KRl ) 

IF  ( ABS ( REAL ( GAMMAF ) ) - GE . ABS ( AIMAG ( GAMMAF ) ) ) 

$GAMMAF=KR1 

IF  ( ABS ( REAL (GAMMAP) ) - GE . ABS ( AIMAG ( GAMMAR ) ) ) 

$GAMMAR=KF1 
GFl(M)=KFl*1.0E+04 
GRl ( M ) =-KRl* 1 . 0E+04 

CALL  DSPRSN (F,U1,U2,E2,D2,SC2 , EPS , GAMMAR , GAMMAF , KR2 , KF2 ) 

KF2=GAMMAF 

KR2=GAMMAR 

GF2 ( M ) =GAMMAF* 1 . 0E+0  4 
GR2 ( M ) =-GAMMAR* 1 . 0E  +  0  4 
CONTINUE 
DO  11  M=1 , NSAMP 
F=*FRQ  ( M ) 

FQ=F* 1 . 0E6 
W=FQ 

ZF1=JJ*UM*W/GF1 ( M) 

ZR1=-JJ*UM*W/GR1 ( M ) 

ZF2=JJ*UM*VJ/GF2  (M) 

ZR2=-JJ*UM*W/GR2 ( M ) 


RM=  (  RMF{  M )  +RflR(  M  )  ) 

XM»(XMF(M)+XMR(M) ) 

ALPH»6.32E-7*SQRT(FQ)/(B0*1.0E-06*ZC)+0.5*RM/ZC 

BETA-3. 6361E-8*FQ+0.5*XM/ZC 

AR1»0.5*RMF(M)/ZC 

AR2-0 . 5*RMR( M ) /ZC 

GAMMA-CMPLX ( ALPH , BETA ) 

ZIN»ZC*CTANH( GAMMA* RLGH*1 .0E-03 ) 

RIN-REAL ( Z IN ) 

RMl-RIN* ARl/ALPH 

RM2-RIN*AR2/ALPH 

VF1-CEXP ( -GF 1 ( M) *B0*1 .0E-04 ) 

VRl-CEXP(-GRl(M)*B0*1.0E-04 ) 

TZl- ( ZFl+ZRl ) 

TZ2-ZF2+ZR2 

C  FIND  T  MATRIX  FOR  METALIZED  SECTION 
AM- ( ZFl*VFl+ZRl*VRl ) /TZ 1 
BM- ( ZFl*ZRl* ( VF1-VR1 ) ) /TZ 1 
CM*( VFl-VRl ) /TZl 
DM- ( ZFl*VRl+ZRl*VFl ) /TZ I 
C  FIND  S-PARAMETERS 

30  CALL  SE(ZIN,RS,ZF2,ZR2,RM2,RMl,SSB) 

CALL  SW(ZF2,ZR2,AM,BM,CM,DM,SS8) 

CALL  PATH (ZF2,ZR2,TZ2,GF2(M) ,A0,GR2(M) ,PP) 

CALL  PATH(ZF2,ZR2,TZ2,GF2(M) ,S1(1) ,GR2(M) ,PPP) 

CALL  CASCAD ( YSB , PPP , YSB , ARR1 ,  F ) 

CALL  CASCAD (YSB, PPP, YSB, ARR2,F) 

33  CALL  FQFP ( ARRI , PP , ARR2 , RS , ZR2 , Z F2 , PLOSS , M ) 

WAV ( M) *2 . *PI/AIMAG ( GF2 ( M ) ) 

11  CONTINUE 

WRITE (6, 87 ) 

87  FORMAT ( ' 1 ' ,10X, ' FRQ , WAV , PLOSS ,TDY '// ) 

DO  88  1  =  1 r NSAMP 

88  WRITE (  6  ,  *  )  I , FRQ ( I ) , WAV { I ) ,PLOSS<I) 

P-0 

CALL  PLOT80 ( FRQ , PLOSS , WAV ,-80 . ,0 . , NSAMP) 

66  FORMAT( 2 5X, 'MATERIAL  PARAMETERS'//'  YIG  THICKNESS 
$= '  F5 . 1 

3'  UM'/' "AIR-GAP  = ' , F5 . 1 ,  '  UM'/'  ALUMINA  *',F7.1,'  UM'/ 
1'  MAGNETIC  BIAS  FIELD  =',F7.2,'  OERSTEDS'/’  MAGNE 
$TIC  SATURATION' 

2 , F7 . 1 , '  OERSTEDS'/'  MINIMUM  LINEWIDTH  =',F5.2,' 
$OERSTEDS ' / ) 

67  FORMAT ( 25X, 'TRANSDUCER  PARAMETERS’//'  THICKNESS 
$-' ,F5.2, 'UM'/ 

3’  WIDTH  * '  , F5 . 1 ,  '  UM'/'  LENGTH  =’,F4.1,'  MM'/' 


n  n 


SSEPERATION  *' 

4 , F7 . 1 ,  '  UM'/) 

68  FORMAT { '  INPUT  TRANSDUCER 15 , ’  BARS;  SPACING 
$  IN  UM: '/15F8.1/) 

69  FORMATC  OUTPUT  TRANSDUCER 15  , '  BARS;  SPACING 
SIN  UM: * /15F8 . 1/) 

61  FORMAT (//25X, 'A  SAMPLE  OF  VALUES'//1  NO.', 15/' 

SFREQUENCY 

1F8.2/'  FREE-YIG  GAMMA  *',2E16.5/*  METALLIZED 
SGAMMA  =' , 2E16 . 5/ 

2'  RADIATION  RESISTANCE  =',E16.5/'  RADIATION 
SREACTANCE  = ' 

3 ,E16 . 5/'  INPUT  IMPEDANCE  »',2F14.5/'  LUMPED 
$RAD -  RES-  *' ,F9-4/) 

78  FORMAT (//3(1X,2F14.5,3X) ) 

STOP 

END 


SUBROUTINE  DSPRSN ( F , U1 , U2 , E , D , S , EPS , GAMMAF , 
$GAMMAR,GF,GR) 

DIMENSION  D ( 5 ) ,  £  (  5  ) ,  S  (  5 ) ,  K  (  5  ) , KP ( 5 ) 

COMPLEX* 8  J,U1,U2,K,G,K12,K32,K132,G34,G54,C1., 
$C2,CX2,CX4,CT3,CT1,K132P,G34P,G45P,C1P,C2P,GR, 

1CT5,TANH13,FG,GMA,DIF,GN,FGMA,KP,K12P,K32P, 

2 , CX2P , CX4P , CT3P , CTlP , CT5  P , TNHl 3P , GAMMAF , GAMMAR 
SFGDG , FGP ,W0 ,GF 

DATA  E0,U0f PYf GYRO/8 .84 19E-1 8 ,1.2566 37E-1 2, 
$3.141592654,2.8/ 

J “CMPLX ( 0 . 0 , 1 . 0  ) 

W*2 -0*PY*F*1 -0E+06 
GMA=GF 

4  DO  5  M*l,700 

CALL  FUNCTN (U0,E0,W,D,E,S,J,PY, GMA ,Ul,U2,K,K12, 
$K32,K132,G34,G54,C1 
1 ,C2 ,CX2 ,CX4 ,CT3 ,CT1 ,CT5 ,TANH13  ) 

FGMA*CX4  * (G54*C2+CT5 ) *(G34*C1+TANH13)-(G54*C1- 
$CT5 ) * ( G34*C2-TANH13 ) 

IF ( CABS ( FGMA ) .LE.EPS)  GO  TO  13 
GN=GMA 

GMA*CSQRT( ( J*CLOG( ( G54 *C1-CT5 ) * ( G34 *C2-TANH 1 3 ) / 
$( (G54*C2+CT5)*(G34* 

1C1+TANH13) ) )/( 2.0*D(4) ) ) ** 2-W*U0* ( W* EO *E ( 4 ) -J* 
$S ( 4 ) ) *(U1**2-U2**2)/ 

2Ul ) 

DIF*GN-GMA 

IF ( CABS (DIF) . LT. ( 1 . 0E-05 *CABS ( GN ) ) )  GO  TO  13 
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IF(M.LE.IO)  GO  TO  5 
GMA* ( GMA+GN ) / 2  -  0 
CONTINUE 

WRITE ( 6 , 10  3 ) M , GN , FGMA ,  DIF , GMA 

FORMAT(10X, 'THERE  IS  NO  CONVERGENCE  IN  ',13,' 

$  ITERATIONS :  FINAL  VAL 

SUES  WERE: '/2  5X, 'GN*' , 2 ( El 5 . 8 , IX ) , 5X , • FGMA*' , 
$2(E15-8,1X)/25X, 'DIF*' 

$ , 2 ( E15 . 8 , IX ) ,5X, 'GMA*' ,2(E15-8,1X>) 

STOP 

GAMMAF=GMA 

WAVE*  2.0*  PY/C ABS ( GAMMAF ) 

GMA*GR 

DO  19  M=1 , 700 

CALL  FUNCTN (U0,E0,W,D,E,S,J,PY, GMA ,U1,U2,K,K12, 
$K32,K132,G34,G54,C1 
1 , C2 , CX2 , CX4 , CT3 ,CT1 , CT5  ,TANH1 3 ) 

FGMA*CX4  * ( G54*C1+CT5 ) * ( G34*C2+TANHl 3 )-(G54*C2- 
$CT5)*(G34*C1-TANH13) 

IF(CABS( FGMA) .LE. EPS)  GO  TO  27 
GN*GMA 

GMA*CSQRT ( ( J*CLOG( (G54 *C2-CT5 ) * { G34*C1-TANH13 ) / 

$( (G54*C1+CT5)*(G34* 

1C2+TANH13 ) ) )/(2.0*D(4) ) ) **2-W*U0 * ( W*E0 *E ( 4 ) -J* 
$S(4))*(Ul**2-U2**2)/ 

2U1) 

DIF*GN-GMA 

IF( CABS ( DIF ) . LT . ( 1 . OE-O 5 *CABS ( GN ) ) )  GO  TO  27 
IF(M.LE.IO)  GO  TO  19 
GMA* ( GMA+GN )/2 .0 
CONTINUE 

WRITE ( 6 , 103 )M,GN , FGMA, DIF, GMA 
STOP 

GAMMAR=GMA 

WAVR*  2.0*  PY/CABS ( GAMMAR ) 

RETURN 

END 

SUBROUTINE  FUNCTN ( U0 ,E0 ,W,D,E,S,J,PI ,GAM,Ul ,U2 , 
$K,K12,K32,K132,G34, 

1G54 ,C1,C2,CX2,CX4 ,CT3,CT1 ,CT5,TANH13 ) 

DIMENSION  D(5) ,E(5) ,S(5) ,K(5) 

COMPLEX *8  X,CSINH,CCOSH,CTANH ,CCOTH ,CCSCH ,CSECH , 
$J,GAM,Ul,U2,K,K12 

1  ,  K32 , K 1 32 ,G34 , G54 , Cl , C2 , CX2 , CX4 , CT3 , CT1 , CT5 , TANHl 3 
CSINH(X)=0.5*(CEXP(X)-CEXP(-X) ) 
CCOSH(X)=0.5*(CEXP(X)+CEXP(-X) ) 

CTANH ( X ) =CS INH ( X ) /CCOSH ( X ) 


CCCTH (X)=CCOSH(X)/CSINH(X) 
CCSCH(X)=2.0/(CEXP(X)-CEXP(-X)  ) 
CSECH(X)=2.0/(CEXP(X)+CEXP(-X) ) 
K(2)=-J*CSQRT(GAM**2+W*U0*(W*E0*E(2)-J*S(2) ) ) 

K  {  4  )  *- J*CSQRT  ( GAM**2+W*U0  *  ( W*  EQ  *E (  4  )  -  J*S ( 4  )  )*( 
$U1**2-U2**2)/U1) 

DO  1  1*1, 5, 2 

K ( I ) =- J*CSQRT ( GAM**2+W*U0  * ( W*EQ  *E ( I ) - J*S ( I) ) ) 
K12*K(1)/K(2) 

K32=K( 3)/K(2) 

K132*K12*K32 

G34=-K  (  3  )  / ( GAM* * 2+W*U0  *U1*  ( W*E0*E  (  4  )  - J*3(  4  )  )  ) 
G54*-K(5)/( GAM**  2+W*U0  *U1 * (W*E0*E(4)-J*S(4) ) ) 
Cl*Ul*K( 4 )+J*GAM*U2 
C2=U1*K( 4 )-J*GAM*U2 
CX2=CEXP(2.0*K(2)*D(2)  ) 

CX4-CEXP(2.0*K(4)*D(4) ) 

CT3=CTANH(K(3)*D( 3)  ) 

IF(D(1) . EQ.1.0E+06)  GO  TO  2 
IF ( D ( 5 ) . EQ.1.0E+06)  GO  TO  5 
CTl=CTANH ( K ( 1 ) *D ( 1  )  ) 

CT5  =  CTAWH ( K ( 5 ) *D( 5  )  ) 

GO  TO  6 

IF(D{5) .EQ.1.0E+06)  GO  TO  3 
GO  TO  4 

CT1=CMPLX(1. 0,0.0) 

CT5*CTl 
GO  TO  6 

CT1=CMPLX(1. 0,0.0) 

CT5=CTANH(K( 5 ) *D( 5) ) 

GO  TO  6 

CT1=CTANH ( K ( 1 ) *D ( 1 ) ) 

CT5  =  CMPLX( 1 .0 ,0 .0 ) 

TANH 1 3* ( ( 1 -0-CX2) *(CT1*CT3+K132 )- ( 1 . 0+CX2 ) * ( K32* 
$CT1+K12*CT3) )/( (1. 

10-CX2)*(CTl+K132*CT3)-( 1.0+CX2 )* (K12+K32*CT1*CT3) ) 
RETURN 
END 

SUBROUTINE  PLOT80  (X , Y  ,  Z , Y1  ,Y2 ,M) 

DIMENSION  X ( M ) ,Y(M) ,Z(M) ,G ( 80 ) ,H(9) 

DATA  BL  ,  ST  ,  DA  ,  VT  ,Rl,R2/lH  II  L  G  '  / 

N=80 

DY* ( Y2-Y1 ) /8 

H ( 1 ) *Y1 

DO  12  K=1 ,8 

H ( K+l ) *H ( K ) +DY 

WRITE (6,29) ( H ( L ) ,L=1,9) 


DO  17  1=1,  M 

NY=IFIX(N*(Y( I)-Y1)/(Y2-Y1 )+0  -  5) 

DI=FLOAT ( I ) /5 - -FLOAT ( 1/5 ) 

DO  20  J=1,N 

DJ  =  FLOAT( J)/10 .-FLOAT ( J/10  ) 

IF ( DI - EQ . 0  - 0 ) GO  TO  25 
G ( J ) =BL 
GO  TO  20 
G ( J ) =DA 

IF ( DJ . EQ . 0 . 0 ) G ( J ) =VT 
CONTINUE 

IF ( NY . GT . N ) GO  TO  15 
IF(NY-LT.l ) GO  TO  16 
G ( NY ) =ST 
GO  TO  17 
G ( N ) =R2 
GO  TO  17 
G( 1 )*R1 

WRITE ( 6,30 )X( I) ,VT,(G(L) ,L»1,N) ,Y(I) ,Z(I) 

FORMAT ( 1X,F7.1 ,2X,A1 ,80A1 ,2F12 .4 ) 

FORMAT( ' 1 ' ,20X, 'PLOT  OF  Y  VERSUS  X  AND  Z'///3X, 
$9 ( FI  0 . 2  )  / 

1  10X  ,  '  1|  *  ,8(9( , '1|  '  )  ) 

RETURN 

END 

SUBROUTINE  TO  COMPUTE  FREQUENCY  RESPONSE 

SUBROUTINE  FQRP ( Yl , PP , Y2 , RS , ZR1 , Z FI , Z IL , M ) 
COMPLEX*8  Yl( 3,3) ,PP( 3,3) ,Y2( 3,3) ,ZR1 ,ZF1,W1,W2, 
$W3,W4,W5,W6 

COMPLEX*8  A , B , C , D , E , NU M , VL S 
DIMENSION  Z IL ( 152) 

GS= 1 .  /RS 


Wl=Yl ( 2 , 2 ) -PP (1,1) 

W2=Y1( l,l)+1.0/ZRl 
W3=Y 1(3,3) +GS 
W4=-PP ( 2,2) +Y2 (1,1) 

W5=Y2 (2,2) +1 .0/ZF1 
W6=Y2 (3,3) +  GS 

NUM=( (W2*Yl(2,3)-Yl(2,l)*Yl(l,3))*/Y2(2,l)*Y2(3,2 
$ ) -Y2 (3,1) *W5 ) ) 

A=W2*Wl*W3+Yl(l,2)*Yl(2,3)*Yl(3,l)+Yl(l,3)*Yl(2,l 
$ ) * Yl ( 3,2) 

B  =  Y  1  (1,3)*Y1(3,1) *W1  +  Y 1 ( 1 ,2 ) *Yl ( 2 , 1 ) *W3+Yl ( 2 , 3  )  * 

$ Yl ( 3,2) *W2 

C=W4*W5*W6+Y2( 1,2)*Y2(2,3)*Y2(3,1)+Y2(1,3)*Y2(2,1 


-  .  < 
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S ) *Y2 ( 3 , 2  ) 

D=Y2(1,3)*Y2(3,1)*W5+Y2(1,2)*Y2(2,1)*W6+Y2(2,3)* 
$Y2( 3,2)*W4 

E=PP ( 1 , 2 ) *PP ( 2 , 1 ) * ( Y1 ( 1 , 3 ) *Y1 ( 3 , 1 ) -W3*W2 ) * ( W5*W6 
$  —  Y  2 ( 2 , 3 ) * Y2 ( 3,2) ) 

VLS=-GS*PP  (  2 , 1 )  *NUM  ' (  A-B  )  *(  C-D)  +E  ) 

Z IL ( M ) =2  0 -*  A LOG 10 ( CABS ( VLS ) ) 

RETURN 

END 

SUBROUTINE  TO  CONVERT  S-PARAtlETERS  TO  Y-PAPAMETERS 
SUBROUTINE  STOY ( S , Z F 2 , Z R2 , RS , Y ) 

COMPLEX *8  S( 3 , 3 ) ,S1,S2,S3 ,A11 ,A12,A13,A21 ,A22  ,A23 
COMPLEX *8  A31 ,A32 , A  33, DET1 , DETIP , DETIN, DET2,DET2N 
$ , DET2P , DET  3 

COMPLEX* 8  DET3P , DET3N ,DEN1 , DEN2 , DEN3  ,  Y ( 3,3)  ,ZF2, 
$ZR2,ZR,ZF,Zl,Z2 
ZR=CSQRT( ZR2 ) 

ZF=CSQRT( ZF2 ) 

Zl=ZR/ZF 
Z2=ZF/ZR 
SI=S( 1 , 1 )+Z2 
S2  =  S( 2,2  )+Zl 
S3=S ( 3 , 3 ) +1  -  0 

DET1P=S(1,1)*S2*S3+S(1,2)*S(2,3)*S(3,1)+S(1,3)* 

$S(2,1)*S(3,2) 

DET1N-S3*S(2,1)*S(1,2)+S2*S(3,1)*S(1,3)+S(1,1)* 

$S(3,2)*S(2,3) 

DET1=DET1P-DET1N 

DET2P=S1*S(2,2)*S3+S( 1,2) *S(2,3)*S( 3,1 ) +S ( 1 , 3 ) * 

$S ( 2 , 1 ) *S ( 3 , 2 ) 

DET2N=S3*S(2,1)*S(1,2)+S(3,1)*S(2,2)*S(1,3)+ 

$S(3,2)*S(2,3)*S1 

DET2=DET2P-DET2N 

DET3P=S1*S2*S( 3,3)+S(l,2)*S(2,3)*S(3,l)+S(l,3)* 
$S( 2,1)*S( 3,2) 

DET3N=S(1,2)*S(2,1)*S(3,3)+S(2,3)*S(3,2)*S1+ 
$S(3,1)*S(1,3)*S2 
DET3=DET3P-DET3N 
A11=(S2*S3-S(2,3)*S( 3,2) )/DETl 
A21-(S( 2,3) *S( 3,1 )-S3*S( 2,1 ) J/DETl 
A31=(S(2,l)*S(3,2)-S(3,l)*S2)/DETl 
A12=(S( 1 , 3)*S( 3,2)-S3*S( 1,2) )/DET2 
A22=(S1*S3-S( 3,1) *S( 1 ,3) )/DET2 
A32=(S( 1  ,2 )*S( 3,  1  )-Sl*3( 3,2) )/DET2 
A13=(3{  1  ,2 )*S( 2,3)-S2*S( 1 ,3) )/DET3 
A23=(S(1,3)*S(2,1;-S(2,3)*S1 )/DET3 
A33=(S1*S2-S( 1 , 2 ) *  S ( 2,1) )/DET3 


RSS=SQRT(RS) 

DEN1=ZF*A11+ZR 
DEN2=ZR*A22+ZF 
DEN3  =  RSS*( A33  +  1 .0  ) 

Y(1 ,l)=(All/ZF-l -0/ZR)/DENl 

Y ( 2 / 1 ) =  (A21* ( 1 -0/ZR+Zl/ZF) )/DENl 

Y(3,1)=2-0*A31/(RSS*DEN1) 

Y( 1 , 2 ) =Al2* ( l-0/ZF+Z2/ZR)/DEN2 

Y  (  2 , 2  )  -  (A22/ZR-1.0/ZF)/DEN2 
Y( 3/2 )=A32*2 -  0/ ( RSS*DEN2 ) 

Y( 1 ,3 ) = Al 3  * ( 1-0/ZF+Z2/ZR) /DEN  3 

Y  (  2 , 3  )  =  A23*(Z1/ZF+1.0/ZR)/DEN3 
Y ( 3 , 3)= (A33-1 -0 )/(RSS*DEN3) 

RETURN 

END 

SUBROUTINE  TO  CALCULATE  THE  WAVE  PORT  PARAMETERS 

SUBROUTINE  SW( ZF1 , ZR1 , A, B , C , D , SS ) 

COMPLEX* 8  ZF1,ZR1,A,B,C,D,SS(3,3) 

COMPLEX*8  SilNUM, SllDEN ,  DETN  , S21NUM, S 21 DEN 
COMPLEX* 8  S22NUM,S22DEN,S12NUM,S12DEN,JJ 
JJ =CM PL X(0. 0,1-0) 

S 1 1NUM=-C*  Z  Fl-D+A+B/ZF 1 

SllDEN=-A*ZRl-D*ZFl+B+C*ZFl*ZRl 

SS ( 1 , 1 ) =SllNUM/SllDEN*ZRl 

DETN=B*C-A*D 

S21NUM=DETN*( ZF1+ZR1 ) 

S21DEN^ZF1*(C*ZR1-D)-(A*ZR1-B) 

SS( 2 ,1 )=S21NUM/S21DEN 

S22DEN=C*ZR1*ZF1-A*ZR1-D*ZF1+B 

S22NUM=D+B/ZR1-C*ZR1-A 

SS( 2 ,2 )=S22NUM/S22DEN*ZF1 

S 1 2NUM=- ( Z  F1  +  ZR1 ) 

S12DEN=ZR1*( C*ZFl-A)-( ZF1*D-B) 

SS ( 1 , 2 ) =S 1 2NUM/S 1 2DEN 

RETURN 

END 

SUBROUTINE  TO  CALCULATE  THE  THIRD  PORT  PARAMETERS 

SUBROUTINE  SE ( Z IN , RS , Z F2 , ZR2 , RM2 , RMl , SS ) 

COMPLEX*8  ZIN,ZF2,ZR2,SS( 3,3) ,JJ , A , B , C 
J J=CMPLX ( 0 . 0 , 1 . 0 ) 

SS( 3,3)*(ZIN-RS)/(ZIN+RS) 

THETA 1=A I MAG ( RM1/Z IN ) /REAL ( RMl/Z IN ) 

THETA2=AIMAG( RM2/Z IN ) /REAL ( RM2/ZIN) 

TA1=ATAN ( THETA  1 ) 

TA2=ATAN ( THETA 2 ) 


U  U  U  U  m 


-W 


C=2 .0*ZIN/( ZIN+RS) 

A= ( RS/ZF2 ) 

B= ( RS/ZR2 ) 

SS(1, 3)«C*CEXP( JJ*TA2)*CABS( ZR2/ZIN) *SQRT(RM2/ 

$ ABS ( REAL ( ZR2 )  )  ) 

SS( 2,3)=-C*CEXP( JJ*TA1)*CABS( ZF2/Z IN ) *SQRT ( RM1/ 

$ABS( REAL( ZF2 ) ) ) 

SS ( 3 , 1 ) -  A*SS (2,3) 

SS(3,2)=B*SS{1,3) 

RETURN 

END 

SUBROUTINE  TO  CASCADE  ZMICROSTRIPS 
G  IS  THE  FIRST  BOX  AND  H  IS  THE  SECOND 

SUBROUTINE  CASCAD ( A , B , C , Y , F ) 

COMPLEX* 8  ABC , CAB , AB , BA , Y ( 3,3),B(3,3),C(3,3), 

$A( 3 , 3 ) ,  ZM , BZMA 

BZMA=B( 1 , 1 ) -A( 2 , 2 ) 

BA  =  B ( 1 , 2 ) *A ( 1 , 2 ) 

AB*B ( 2 , 1 ) *A( 2 , 1 ) 

ZM-C ( 1 , 1 ) -B ( 2 , 2 ) +B ( 2 , 1 ) *B ( 1 , 2 ) /BZMA 
ABC= (  B  (  2 , 1 )*A{2,3) /BZMA-C ( 1 , 3 ) ) /ZM 
CAB  = ( C ( 3 , 1 ) -A ( 3 , 2 ) *B ( 1 , 2 ) /BZMA ) /ZM 

Y(1,1)=A(1,1) +A (2,1)*A(1,2) /BZMA-BA* AB/ ( BZMA*BZMA*  ZM ) 
Y(1,2)=B(1,2)*C(1,2)*A(1,2)/(BZMA*ZM) 
Y(1,3)=A(1,3)+A(1,2)*A(2,3) /BZMA-BA/BZMA* ABC 
Y ( 2 , 1 ) =  B(2fl)*C(2,l)*A(2,l)/( BZMA*  ZM ) 
Y(2,2)*(C(2,2)-C(2,1)*C(1,2)/ZM) 

Y(2r3)=(C(2,3)+C(2,l) *ABC ) 

Y( 3,1 )=A( 3,1 )+A( 3,2)*A( 2 , 1 ) /BZMA+AB*CAB/BZMA 
Y(3,2)=C(3,2)-C(1,2)*CAB 

Y( 3 ,3)=C( 3 ,3)+A( 3 ,3)+A( 3 ,2) *A( 2 , 3 ) /BZMA+ZM*CAB*ABC 
WRITE ( 6 , 34 ) 

4  FORMAT ( '  MIN') 

C  WRITE ( 6 , * ) F , Y 

RETURN 
END 

C  SUBROUTINE  TO  FIND  Y-PARAMETERS  OF  LAND  REGION 
SUBROUTINE  PATH ( Z F 2 , Z R2 , TZ 2 , GF 2 , AO , GR2 , PP ) 

COMPLEX*3  ZF2,ZR2,TZ2,GF2,GR2,PP( 3,3) ,WF,WR,AP, 

$BP , CP , DP 

38  WF=CEXP(-GF2*A0*1 -0E-04 ) 

WR=CEXP ( -GR2* AO  *  1 . 0E-04  ) 

AP  =  (  Z  F2  *WF  +  Z R2*V,JR )  /TZ  2 
BP  =  (ZF2*ZR2*( WF-WR ) ) /TZ 2 
CP= ( WF-WR ) /TZ  2 
DP=( ZF2*WR+ZR2*WF)/TZ2 
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PP(1,1)=  AP/BP 

PP(1, 2)— 1-0/BP 

PP(  2,1)  — (AP*DP-BP*CP)/BP 

PP (2,2)=  DP/BP 

RETURN 

END 
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3.0  Variable  Time  Delays  Using  Cascaded  Linearly  Dispersive  Delay  Lines 

3.1  Linearly  Dispersive  Time-Delay  Control  Of 
Magnetostatic  Surface  Wave  By  Variable 
Ground  Plane  Spacing 

3.1.1  Introduction 

The  propagation  characteristics  of  magnetostatic  waves  in 
expitaxial  yttrium  garnet  (YIG)  films  have  been  extensively 
investigated  both  experimentally  and  theoretically  for  a 
number  of  years  (1-4).  These  studies  have  shown  that 
the  delay  versus  frequency  characteristics  of  MSSW  delay 
lines  of  a  YIG  film  spaced  from  a  ground  plane  are  function 
of  the  saturation  magnetization,  bias  field,  YIG  film 
thickness,  and  ground  plane  spacing.  A  study  of  how 
different  parameters  affect  the  delay  characteristics  of 
surface  waves  indicated  that  dispersion  control  by  ground 
plane  spacing  is  most  effective  (4).  This  chapter  described 
the  theoretical  and  experimental  results  on  the  theoretical 
simulation  and  experimental  evaluation  of  a  magnetostatic 
surface  wave  variable  ground  plane  spacing  delay  line  with 
linear  group  delay. 

3.1.2  Theoretical  Simulation 

Fig.  3.1  shows  the  variation  of  the  group  delay  with 
frequency  for  MSSW  delay  lines  with  ground  plane  spacing 
(Ti)  and  propagation  length  (Li)  as  adjusting  parameters. 

The  YIG  film  thickness  used  in  these  calculations  is  30  uM. 
Examination  of  this  figure  indicates  that  it  is  possible  to 
achieve  a  linear  variation  of  delay  by  adjusting  the  length 
of  the  short  section  (Li)  of  the  step  ground  plane  as  shown 
in  Fig.  3.2.  A  nine  step  ground  plane  structure  is  selected, 
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Fig-  3»1  Variation  of  group  delay  with  frequency 
for  MSSW  delay  lines  with  ground  plane 
spacing  (Ti)  and  propagation  distance 
(Li)  as  adjusting  parameters. 
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Fig.  3.2  MSSW  Delay  line  wi 
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with  a  eight  steps  ground  plane  in  the  middle  of  this  MSSW 
steped  ground  plane  delay  line.  The  separation  of  this 
eight  steps  ground  to  the  YIG  film  can  be  adjusted  by  a  pair 
of  set  screw  .  The  nineth  step  in  this  MSSW  delay  is  fixed 
by  the  substrate  thickness  which  is  250  uM  with  propagation 
length  can  be  adjusted  by  changing  the  single  bar 
transducer  separation  to  the  edge  of  the  substrate 
photographically.  The  length  of  each  section  is  then 
optimized  to  provide  a  minimun  deviation  from  quadratic 
phase  over  the  linearized  delay  bandwidth.  The  linear  time 
delay  characrter istic  which  is  obtained  by  summing  the 
individual  time  delay  corresponding  to  each  section  of  the 
MSSW  9  stepped  ground  plane  structure,  is  shown  in  Fig.  3.3 
A  linear  regression  fit  is  also  shown  in  this  figure  which 
indicated  a  linear  delay  region  between  2.89  to  3.33  GHz 
with  a  positive  slope  of  26.6  nsec/100  MHz  for  the  nine 
stepped  ground  plane  structure  shown  in  Table  3.1.  The 
calculated  phase  error  in  degrees  from  a  quadratic  phase  of 
this  linear  delay  region  between  2.89  to  3-33  is  shown  in 
Pi9-  3.4.-  The  R.M.S-  deviation  from  quadratic  phase  over 
this  440  MHz  bandwidth  is  8  degrees.  The  corresponding  time 
bandwidth  product  (TBP)  is  near  120  for  this  nine  stepped 
ground  plane  MSSW  theoretical  simulated  linear  delay  line. 
3.1.3  Experimental  Evaluation 

Due  to  difficulty  and  long  waiting  period  for  the 
fabrication  of  the  eight  stepped  centre  ground  plane,  a 
continuously  variable  ground  plane  is  used  to  provide  the 
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3.3  The  summation  of  the  individual  time  delay  section 
in  the  9-step  ground  plane  structure,  with 
dimension  shown  in  table  2.1,  results  a  linear 
time  delay  characteristic  between  2.89  to  3.33 
GHz.  Thickness  of  YIG  is  30  uM. 


Ti(uM) 


Li (cm) 


300 

.18 

250 

.20 

225 

.15 

200 

.15 

175 

.23 

150 

-22 

135 

.20 

115 

.10 

100 

.07 

Table  2.1 


Ground  plane  spacings  and  lengths  for  a  nine 
section  delay  line 


FREQUENCY  (  GHz  ) 


3.4  The  quadratic  phase  error  of  the  linear  delay 
region  of  the  9-step  ground  plane  structure 
The  R.M.S.  deviation  from  quadratic  phase  in  this 
bandwidth  is  8  degrees. 


continuous  variation  of  ground  plane  spacing  to  the  YIG  film 
along  the  propagation  direction  of  the  MSSW  delay  line. 
The  separation  of  this  continuously  variable  ground  plane  to 
the  YIG  film  can  be  adjusted  by  changing  structure  of  the 
ground  plane  and  varying  the  distance  of  this  ground  to  the 
film  by  a  pair  of  set  screw. 

Fig-  3.5  shows  a  sketch  of  the  MSSW  variable  ground  plane 
delay  line.  The  device  is  in  a  "flipped"  configuration  with 
the  YIG  film  flipped  over  the  single  bar  microstrip 
transducer.  A  pair  shorted-circuited  50  uM-wide,  3  mm-long, 
and  3  uM-thick  microstrip  couplers  are  used  for  the  lunching 
and  receiving  of  the  microwave  energy  for  MSSW  propagation 
along  the  YIG  film,  a  25  uM-thick  3  mm-wide  and  15  mm-long 
YIG  film  grown  by  liquid  phase  epitaxy  On  a  500  uM-thick 
gadolinium  gallium  garnet  (GGG)  substrate  was  used  for  the 
delay  measurement.  The  HP  8409B  automatic  network  analyzer 
was  used  to  measure  the  reflection  and  transmission 
parameters  of  the  device  under  test. 

The  group  delay  was  then  obtained  from  the  S12  phase  data  by 
calculation  of  phase  slope.  Since  this  is  a  differentiation 
process,  the  resolution  is  determined  both  by  the  frequency 
aperture  (  f)  and  the  phase  detector  resolution  of  the 
network  analyzer.  By  increasing  the  frequency  aperture  and 
the  corresponding  phase  different,  an  averaging  effect  is 
performed  to  reduce  the  time  delay  ripples.  A  30  MHz 
frequency  aperture  is  chosen  in  order  ot  optimize  the 
calculated  time  dealy  obtained  from  the  automatic  network 
analyzer  phase  measurement.  A  modification  of  the  control 


program  of  the  automatic  network  analyzer  is  performed  to 
computed  the  absolute  phase  change  as  a  function  of 
freguehcy  over  the  frequency  bandwidth  of  interest.  For  a 
linear  delay  bandwidth,  the  corresponding  absolute  phase 
function  is  a  quadratic  function  of  frequency  (Fora 
constant  delay  bandwidth,  the  corresponding  absolute  phase 
fuction  is  a  linear  function  of  frequency).  A  root  mean 
square  fit  is  then  performed  for  this  experimental 
determined  phase  data.  The  deviation  from  quadratic  phase 
defined  as  the  difference  between  the  quadratic  fit  phase 
function  to  the  experimental  phase  data  is  then  used  as  the 
parameters  to  determine  the  deviation  from  linear  delay 
(quadratic  phase)  over  the  bandwidth  of  interest  of  the 
device  under  test. 

By  carefully  adjusting  the  separation  of  the  continuously 
variable  ground  plane  to  the  YIG  film  ,  linear  delay  of 
bandwidth  up  to  500  MHz  can  be  experimental  obtained  with 
phase  error  less  than  20  degrees.  The  slope  of  the  linear 
region  can  also  be  adjusted  by  changing  the  structure  of  the 
continuously  variable  ground  plane.  Fig.  3.6  shows  the 
experimental  delay  measurement  of  one  of  the  device  under 
test.  A  linear  delay  of  500  MHz  bandwidth  (between  2.62  to 
3.12  GHz)  with  a  slope  of  15.3  nsec/100  MHz  was  observed. 
Fig.  3.7  shows  the  quadratic  phase  error  of  this  device  in 
the  linear  delay  region.  The  corresponding  R.M.S.  phase 
error  is  13  degrees  The  insertion  loss  over  the  linear 
region  is  between  15  to  20  dB  (unmatched).  The  separation  of 
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(GHi) 

The  quadratic  phase  error  of  the  continuously 
variable  ground  plane  MSSW  delay  line 
The  R.M.S.  phase  error  is  13  degrees  over  the 
linear  delay  bandwidth. 


this  continuously  variable  ground  plane  to  the  YIG  is 
roughly  determined  by  removing  the  "flipped"  over  YIG  film 
and  measuring  the  profile  of  the  MSSW  delay  line  along  the 
direction  of  propagation.  A  50  uM  air  gap  was  added  to  the 
measuring  data  to  represent  the  air  gap  which  is  unavoidable 
in  the  flipped  configuration.  A  15-step  ground  plane 
structure  ,with  dimensions  given  is  Fig.  3.5  /  was  then  used 
to  approximate  the  continuously  variable  ground  plane 
structure  measured.  The  computing  delay  characteristic  of 
this  15-step  model  is  also  in  Fig.  3i6.  The  difference 
between  this  computing  delay  characteristic  to  the 
exper timental  linear  delay  character ictic  is  quite  large. 
This  large  discrepancy  is  due  to  the  error  associated  with 
the  measurement  of  the  bias  field  strength,  YIG 
th ickness  ,and  separation  of  the  ground  plane  to  the  YIG 
film . 

Applications  of  these  MSSW  linear  delay  devices  include 
pulse  compression  system  (5,6),  compressive  receivers  and 
variable  time  delay  devices  utilizing  a  cascade 
configuration  with  a  linearized  backward  volume  wave  device 
(7,8)  . 
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3.2  Electronically  variable  Time  Delays  Using  Cascaded 
Magnetostatic  Wave  Delay  Lines 

3.2.1  Introduction 

Two  features  which  make  magnetostatic  wave  (MSW)  devices 
attractive  are  (a)  an  operating  frequency  range  above  one 
GHz  and  (b)  electronic  tunability.  Both  of  these  features 
are  required  for  electronically  variable  time  delay 
components  used  in  phased  array  radar  and  communication 
systems.  The  basic  architecture  for  variable  time  delay 
device  based  on  MSW  technology  was  described  originally  by 
Setheres,  Owens  and  Smith  (7).  A  Schematic  of  their  approach 
is  illustrated  in  Figure  3.3.  Briefly,  two  MSW  delay  lines, 
one  biased  to  operate  in  the  backward  volume  wave  (MSBVW) 
mode  and  the  other  biased  to  operate  in  the  surface  wave 
(MSSW)  mode,  are  cascaded  together.  The  dispersion  (i.e., 
the  delay  vs.  frequency)  of  the  MSSW  device  has  a  positive 
slope,  while  the  slope  of  the  MSBVW  dispersion  is  negative- 
If  the  absolute  magnitudes  of  these  two  slopes  are 
identical  ,  and  if  the  dispersions  are  linear,  then  the  net 
dispersion  of  the  cascaded  .device  will  be  zero  over  a  common 
frequency  band.  If  the  frequency  domain  of  one  delay  line  is 
shifted  by  changing  the  strength  of  its  applied  biasing 
field,  then  the  absolute  delay  of  the  cascaded  device  will 
also  be  shifted.  Thus,  different  values  of  constant  delay 
can  be  obtained  by  simply  changing  the  strength  of  the  bias 
by  a  small  amount.  One  impediment  to  the  realization  of  a 
device  of  this  with  large  bandwidth  is  the  fact  that  in 
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general  the  dispersion  of  typical  MSW  devices  is  non-linear. 
In  the  original  experiments  no  attempt  was  made  to  linearize' 
the  individual  dispersions.  Even  so  the  results  were 
generally  encouraging.  In  this  chapter  an  improved  version 
the  the  original  device  in  which  care  has  been  taken  to 
linearize  the  dispersion  for  both  the  surface  and  backward 
modes  is  reported. 

3.2.2  Theory 

From  the  study  of  the  dispersive  characteristics  for  MSSW 
and  MSBVW  delay  line,  it  is  found  that  the  shape  of  the 
dispersion  curve  is  strongly  dependent  on  the  thickness  of 
the  film  and  the  spacing  above  the  ground  plane.  It  has  been 
demonstrated  (9)  that  the  linearity  of  the  dispersion  can 
be  substationally  improved  by  carefully  adjusting  the' 
spacing  between  the  ground  plane  and  the  YIG  film.  Thus,  the 
ground  spacing  technique  was  incorporated  into  the  design  of 
the  delay  line  components  of  the  cascaded  device - 
A  linear  MSBVW  dispersion  curve  can  be  obtained  in  a 
straight  forward  manner  by  using  a  35  uM  thick  layer  of  YIG 
spaced  250  uM  above  the  ground  plane.  This  spacing  is  very 
close  to  the  standard  for  MSW  delay  lines  and  no  special 
fabrication  techniques  are  required,  other  than  some 
polishing  of  the  dielectric  spacer.  It  requires  a  much 
closer  ground  plane  spacing  to  achieve  the  same  results  with 
the  MSSW  mode  and  past  experience  has  shown  that  ground 
current  losses  increase  as  the  spacing  decreases. 

To  circumvent  this  problem  a  variable  ground  plane  structure 


of  the  type  reported  in  section  by  Chang,  Owens  and  Carter 
(10)  was  selected  for  the  MSSW  device.  The  essence  of  this 
approach  is  that  the  YIG  delay  line  is  held  above  a  ground 
plane  whose  distance  from  the  YIG  film  is  a  function  of  the 
path  length.  A  schematic  of  the  technique  is  shown  in  Fig. 
3.5.  The  advantage  of  this  approach  4s  that  in  principle 
optimum  dispersion  control  can  be  achieved  and  the  minimum 
spacing  required  is  greater  than  in  the  uniform  spacing 
case . 

For  both  modes  of  propagation  the  figure  of  merit  is  defined 
in  terms  of  the  phase  deviation  from  the  appropriate 
theoretical  phase  vs.  frequency  curve.  If  the  delay  time  is 
a  linear  function  of  the  frequency  then,  the  phase 
dependence  on  frequency  is  quadratic.  On  the  other  hand,  if 
the  delay  time  is  constant  over  the  frequency  band,  then  the 
phase  is  a  linear  function  of  frequency.  In  principle,  it  is 
possible  to  achieve  close  to  zero  deviation  from  quadratic 
phase  over  a  bandwidth  of  500  MHz  centered  at  3  GHz  using 
the  ground  plane  spacing  tephniques  described  in  chapter  2. 

3.2.3  Experimental  Results 

A  cascaded  device  was  prepared  using  a  MSSW  delay  line  with 
a  variable  ground  plane  and  a  MSBVW  delay  line  with  a  YIG 
film  thickness  35  microns  and  ground  plane  spacing  of  254 
microns.  The  characteristic  of  the  MSBVW  delay  line  is  shown 
in  Fig.  3,9,  The  variable  ground  plane  was  prepared  by 
polishing  a  ground  plane  into  a  smooth  curve.  The  ground 
plane  curve  was  optimized  by  alternating  vs.  frequency 


measurements  with  block  sanding  until  the  desired  dispersion 
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ria  39  The  characteristic  of  the  MSBVW  delay  line 
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curve  was  obtained. 

The  delay  line  material  was  prepared  by  growing  YIG  films  on 
one  inch  diameter  gallium  gadolinium  garnet  (GGG)  wafers 
using  the  liquid  phase  epitaxy  (LPE)  method.  The  YIG 
thickness  for  the  MSSW  delay  line  was  25  microns,  while  the 
thickness  of  the  MSBVW  delay  line  YT'j  was  35  microns.  The 
wafers  were  cut  into  bars  for  delay  line  use,  and  the  ends 
of  these  bars  were  beveled  to  reduce  reflections. 

The  individual  delay  lines  were  experimentally  evaluated 
with  an  automatic  network  analyzer.  Large  laboratory 
electromagnets  were  used  to  supply  the  necessary  d.c.  bias. 
In  the  final  conf iguration  both  the  MSSW  delay  line  and  the 
MSBVW  delay  line  exhibited  R.M.S.  phase  deviation  from 
quadratic  of  about  13  degree  over  approximately  500  MHz 
bandwidth  centered  at  3  GHz.  Experimental  plots  from  the 
network  analyzer  of  delay  time  vs.  frequency  for  both  MSSW 
and  MSBVW  devices  are  given  in  Figure  3.10.  The  slopes  of 
these  curves  are  9.6  nsec/100  MHz  and  -9.7  nsec/100  MHz, 
respectively. 

It  is  found  that  by  adjusting  the  d.c.  bias  field  of  the 
virable  ground  plane  MSSW  delay  line,  the  dispersion 
characteristic  of  the  linear  delay  bandwidth  becomes 
nonlinear.  In  order  the  deviation  from  quadratic  phase  over 
the  frequency  band  of  interest  increases  substantially  by 
changing  the  bias  field.  For  the  backward  volume  wave  delay 
line  although  the  linear  bandwidth  is  changed  by  adjusting 
the  bias  field  but  the  dispersion  characteristic  over  the 
bandwidth  of  interest  remains  linear  with  a  slightly 


adjustment  in  the  slope  Table  3.2  given  the  effect  of 
changing  the  d.c.  bias  of  the  MSBVW  delay  line  on  the  slopes 
of  the  linearized  frequency  band.  Because  of  these 
experimental  results,  the  bias  field  for  the  MSSW  linear 
delay  line  in  fixed  at  the  same  value  whereas  the  bias  field 
of  the  MSBVW  linear  delay  line  can  be  changed  to  provide  the 
adjustment  for  the  time  delay  in  the  cascaded  device. This 
two  MSW  delay  line  is  then  cascaded  together  to  performed 
the  variable  time  delays  experiments. 

Delay  vs.  frequency  curves  experimentally  obtained  from  the 
cascaded  MSSW  and  MSBVW  device  are  shown  in  Pig.  3.11.  To 
obtain  these  curves  the  bias  on  the  MSSW  delay  line  was  held 
constant  at  410  Gauss  while  the  bias  applied  to  the  MSBVW 
delay  line  was  varied  from  540  to  650  Gauss.  The  data 
obtained  are  summerized  in  Table  3.3:.  The  maximum  delay 
differential  obtained  was  42  nsec  with  an  applied  bias 
differential  of  110  Gauss.  The  bandwidth  of  this  variable  is 
between  2.9  to  3.15  GHz  ,  and  the  R.M.S.  phase  deviation 
from  linearity  ranged  from  6.4  degrees  to  12.9  degrees. 
These  values  are  quite  acceptable  for  a  number  of  phased 
array  radar/communicatins  applications  and  represent  a 
significant  advancement  in  the  state  of  the  art  for  MSW 
devices . 

There  is  one  problem  with  the  cascaded  device  as  it  stands 
now.  The  insertion  loss  is  high,  approximately  35  db  across 
the  band.  This  high  value  is  due  in  large  part  to  the  fact 
that  the  device  was  not  tuned,  and  the  loss  should  be 
reduced  consederably  with  proper  matching  of  the 
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Table  3,2  The  effect  on  changing  the  slopes  of  the  MS3VW 
linear  delay  line  by  adjusting  the  bias  field 
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Fig.  3.11  Dispersion  curves  for  ca  scad  delay  line  as  a 
function  of  bias  field 
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Table  '3.3  Experimental  results  for  the  MSW  Variable  Time 
delay  lines 


transducers. 


3,2.4  Summary 

An  electronically  variable  time  delay  device  comorised  of 
two  cascaded  MSW  delay  lines  has  been  fabricated  and 
evaluated,  and  the  operating  characteristics  are 
substationally  better  than  previously  reported  embodiments 
of  this  device.  The  improvements  are  due  mainly  to  the 
linearization  of  the  two  component  dispersions.  With  a 
reduction  in  insertion  loss  and  compact  packaging  a  device 
of  this  type  will  be  ready  for  service  in  phased  array 
systems . 
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4.0 


3-6  GHz  Tunable  MSSW  Resonator 


Magnetostatic  wave  propagation  in  EPI  YIG  films  satisfy  the  equation 
V  x  h  =  0  or  h  =  Vi|i. 

The  wave  can  propagate  in  three  modes  such  as  Surface,  Forward  volume 
wave  and  Backward  volume  wave  (Fig.  4.2) 

The  surface  wave  mode  is  excited  when  the  bias  field  lies  on  the  surface 
of  the  YIG  films  such  that  the  field  direction  is  normal  to  the  wave 
propagation  direction. 

For  Backward  volume  wave,  the  field  is  in  the  direction  of  the  wave 
propagation  and  for  forward  volume  wave,  the  field  is  perpendicular  to 
the  YIG  surface. 


INTRODUCTION 

Magnetostatic  wave  (MSW)  devices  form  the  basis  of  a  new  emergence  of 
microwave  analog  signal  processing.  Among  these  devices,  MSW  resonators 
offer  over  surface  acoustic  wave  (SAW)  resonators  simple  transducers  not 
needing  narrow- 1 ine-width  photolithography  and  broad  frequency  tunability 
by  an  externally  applied  magnetic  field.  Over  yttrium  iron  garnet  (YIG) 
sphere  resonators  which  operate  well  at  microwave  frequencies,  but  which 
require  tedious  and  expensive  fabrication  procedures,  MSW  techniques 
provide  means  of  obtaining  high  Q  resonators  allowing  MSW  oscillator  and 
complex  filter  functions. 

MSSW  are  predominantly  magnetic,  slow,  dispersive  electromagnetic  waves 
(with  dispersion  shown  above)  that  propagate  with  the  field  maxima  at  the 
ferrite  surfaces.  For  the  forward  propagation  mode,  if  the  wave  energy 
is  confined  to  the  top  surface  of  the  YIG  then  for  the  reverse 
propagation  mode,  the  energy  is  confined  to  the  bottom  surface  of  the 
YIG. 

The  basic  element  of  a  MSSW  Resonator  is  a  periodic  grating  on  the  YIG 
surface  which  reflects  MSSW  at  specific  wavelength  determined  by  the 
period  of  the  grating. 

A  schematic  of  a  two-port  MSSW  resonator  is  shown  in  figure  4.6  The 
two-port  resonator  consists  of  a  pair  of  distributed  mirrors  separated  by 
a  propagating  spacing.  The  input  and  output  loop  transducers  are  usually 
placed  within  the  spacing  between  the  two  arrays  with  the  structure  (a) 
on  top  of  the  structure  (b)  in  flipped  configuration.  The  structure 
provides  an  easy  isolation  between  the  input  and  the  output  port. 
F.fficient  distributed  reflectors  are  formed  by  several  types  of 
discontinuities  (1)  chemically  or  ion  beam  milled  etched  groove  arrays  or 
(2)  metal  dot  arrays.  In  this  work,  the  grooves  are  plasma  etched,  and 
theresonators  were  designed  to  operate  at  A  =  300  pm. 

In  this  work,  we  have  attempted  on  the  improvement  of  the  following 
characteristics  of  a  MSSW  resonator. 


;  FIGURE  4-.li  Dispersion  Relations  of  Magnetostatic 

Waves  in  YIG  Plates 


Figure  4.3 


Figure  4.4 


Figure  4. 


Dispersion  curves  for  the  case  of  magnetostatic  wave 
propagation  along  the  biasing  field.  The  curves  designated 
by(a) , (b) ,  (c) ,  etc.,  correspond  to  n=0,l,2,  etc.  The  slope  of 
U)  vs  k  curves  is  negative  throughout  the  frequency  range  for  all 
mode  orders.  Thus  the  modes  are  backward  waves  for  which  the 
directions  of  energy  flow  and  phase  propagation  are  mutually 
opposite  (after  Damon  and  Eshback,  1961). 


<■>)  (b) 

Magnetostatic  dispersion  curves  for  (a)  free  ferrimagnetic 
slab,  and  (b)  metal-backed  ferrimagnetic  slab.  In  (b)  S=±l 
corresponds  to  propagation  along  the  ±y-direction.  The  introduction 
of  the  metal  plate  thus  removes  the  degeneracy  and  renders  the 
propagation  characteristics  nonreciprocal  (after  Seshadri,  (1970). 


Magnetostatic  mode  spectrum  for  a  normally  magnetized  slab. 
The  phase  velocity  is  smaller  for  the  higher-order  modes 
(after  Damon  and  van  de  Vaart,  1965). 
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Figure  4.6  Schematic  sketch  of  a  resonator  and  transduce 
combination 


16  5 


Figure  4.7  Resonator  (#  K2)  response  at  6  GHz  with  the  matched  1  mm 
transducer.  The  bottom  trace  correspond  to  The  dip 

near  resonance  shows  matching.  The  insertion  loss  at 
resonance  is  14  dB.  It  can  be  seen  that  the  insertion 
loss  can  be  further  improved  if  the  resonance  is  placed 
at  the  center  (top)  of  the  transducer  passband.  This  will 
increase  the  out  of  resonance  rejection. 


Figure  4,8  Resonator  response  (/>  K2)  showing  13  dB  insertion  loss  at  6 
GHz.  1  mm  matched  transducer  is  used. 
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1)  Smooth  tunability  of  resonance  over  3-6  GHz 

2)  High  loaded  quality  factor  Q  without  the  necessity  of  insertion  of 
any  dielectric  decoupler  between  the  YIG  and  the  transducer  as  it 
increases  the  insertion  loss. 

3)  Insertion  loss  lower  than  those  obtained  so  far  in  literature 
(CASTERA  et  al) 

4)  High  off  resonance  rejection  and  single  mode  operation 
EXPERIMENTAL  RESULTS 

Our  works  show  that  it  is  possible  to  meet  all  the  above  four  criteria  in 
1  mm  wide  resonators  built  with  EPI  YIG  (YIG  thickness  d  14  pm,  etch 
depth  h/d  =  2.5%  #  bars  =  65,  array  separation  2.4  mm,  L  =  75  pm,  GGG 
thickness  =  500  pm)  and  using  matched  loop  transducer  with  1  mm  long, 
microstrip  built  on  gold  plated  quartz  substrates  of  thickness  250  pm  (50 
pm  strip  width) . 

The  measured  loaded  Q  lies  between  600-2000  over  3-6  GHz  with  insertion 
loss  8-14  dB  (13  dB  IL  at  6  GHz  has  also  been  obtained  for  one  resonator 
Fig.  47  The  off  resonance  rejection  is  10-l4dB,  higher  at  high  fre¬ 
quencies.  The  saturation  power  is  4  dBm  at  3  GHz  and  increases  with 
frequency.  The  resonators  are  tunable  over  3-6  GHz  (see  the  attached 
figures) . 

As  is  apparent  from  the  trace  (Figure  4 . 9 ) that  sidelobe  rejection  might 
possibly  be  increased  to  20  dB  or  better  by  carefully  adjusting  the 
separation  of  the  two  bars  of  a  loop  transducer.  This  will  change  the 
wavelength  of  the  launched  wave,  which  can  shift  the  resonance  peak 
exactly  on  top  of  the  transducer  response.  This  later  improvement  is 
being  attempted. 

3  mm  wide  and  1.5  mm  wide  resonators  were  also  fabricated  along  with  3  mm 
long  and  1.5  mm  long  loop  transducer. 

Loaded  quality  factor  Q  in  the  range  800-2800  over  3-6  GHz  was  obtained 
with  3  mm  resonator.  (The  separation  between  the  two  arrays  ©1.5  mm.) 

The  corresponding  insertion  loss  was  15-34  dB.  The  off  resonance 
rejection  is  found  to  be  only  6  dB.  However,  the  resonators  in  general 
are  not  usable  beyond  4.5  GHz  because  of  the  appearance  of  the  uniform 
resonance  mode  higher  than  the  MSSW  resonance,  (at  low  frequency  side  of 
the  MSSW  resonance)  Q  in  general  increases  at  higher  frequencies  so  does 
the  insertion  loss.  (A  typical  3  mm  resonator  response  is  shown  in  Fig. 
4.11) 

Using  1.5  mm  wide  MSSW  resonators  and  1.5  mm  long  transducer,  we  found 
that  MSSW  resonance  is  tunable  over  3-6  GHz,  however  the  insertion  loss 
16-29  dB  was  still  considered  to  be  high.  The  corresponding  Q  was  found 
to  be  lower  than  3  mm  resonator  (namely  650-2400). 

In  order  to  reduce  the  insertion  loss  further  1  mm  wide  resonators  were 
made  along  with  1  mm  gold  transducers  on  quartz  substrates,  though  this 
reduced  Q  slightly  (600-2000  over  3-6  GHz),  the  insertion  loss  was 
reduced  appreciably  to  8-20  dB. 

After  it  is  found  that  using  narrow  (1  mm  wide)  resonator  and  narrow 
transducer  strip  (1  mm  long  microstrip)  the  insertion  loss  was  reduced  to 
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8-20  dB,  it  was  decided  to  match  the  transducer  to  50  ohm  input  output 
load.  Sjj  parameters  (both  magnitude  and  phase)  were  determined  at 


resonance  over  3-6  GHz. 


To  measure  S. .  calibration  was  done  over  3-6  „„ 

1 1  GHz 


with  the  corresponding  transducer  short  i.e.,  by  shorting  the  1  mm 
transducer  feeder  line  at  the  input  of  the  loop  microstrips. 

Sjj  parameters  were  plotted  on  Smith  Chart  and  primary  matching  design 

principles  were  followed  as  provided  by  Hewlett  Packard.  (Application 
note  154).  As  Matching  over  an  octave  band  3-6  GHz  is  difficult,  we 
tried  to  optimize  the  matching  through  COMPACT  (computer  optimization  of 
microwave  passive  and  active  circuits).  The  measured  parameters  over 

3-6  GHz  are  shown  in  Smith  Chart  polt  A.  The  circuit  which  optimized  the 
matching  is  shown  below:  (To  be  built  on  250  pm  quartz  substrates 


This  circuit  is  found  to  reduce  the  S ^ ^  parameters  within  4000-5850  mHz 

as  shown  in  Smith  Chart  plot  B  which  shows  the  computer  optimized  S^ 

parameters.  We  were  not  much  concerned  about  optimization  in  the  range 
3-4  GHz,  because  the  insertion  loss  at  resonance  was  within  8-12  dB  even 
with  an  unmatched  1  mm  transducer. 

Microstrip  Matched  Circuit: 

Inductance  and  capacitance  of  the  computer  optimized  circuits  were 
realized  through  the  following  circuit  equations 

1)  High  Impedance  Microstrip  Lines 

Impedance  around  100  ohm  were  used.  (ZoL  =  100Q) 

Required  length  of  the  microstrip  line  for  the  inductance  (1.5  nH) 
was  determined  as  described  below 

2)  Discontinuity  and  inductance  were  utilized  in  the  circuit  as  shown 
in  Fig.  4.12 
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4.9  MSSW  resonator  with  1  mm  unmatched  transducer  at  6  GHz. 

The  resonator  is  designed  to  operate  at  300  micron 
wavelength.  The  insertion  loss  is  20  dB  with  Q  at  resonance 
2000.  The  upper  trace  corresponds  to  S^.  It  is  seen  that 

at  the  resonance  point  there  is  not  a  significant  dip  in  S 
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Figure  4. 10  MSSW  resonator  response  with  1  mm  matched  transducer  at 

6  GHz.  The  insertion  loss  at  resonance  is  \k  dB  with  the 
same  resonator  as  above.  The  bottom  trace  corresponds  to 
Sjj.  It  is  seen  that  there  is  a  significant  dip  in  a 

resonance.  The  at  resonance  is  0.15.  The  resonator 

identification  number  is  K2. 


;  j.sy.i,  ■ 


Figure  4. 11A  typical  3  mm  wide  resonator  response  with  3  mm 
transducer  at  5.5  GHz. 


Resonator  identification  #6 
f  =  5.5  GHz 
Af  =  1  GHz 


Figure  4. 12  Matched  input  transducer  circuit  as  obtained  through 
compact 


The  frequency  in  the  following  equations  were  chosen  to  be  the 
center  of  the  band  (3-6  GHz). 
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£efj  for  inductor  microstrip  AB  (Zol  =  100  Q)  can  be  obtained  from  equa¬ 
tion  (3). 


£ej.fr  for  capacitance  microstrip  can  be  determined  by  first  finding  out 

Z  corresponding  to  the  microstrip  of  width  'w'  =  0.55  mm  from  the 

equations . 
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using  the  corresponding  eg££  as  obtained  above. 

As  applicable  to  the  microstrips  and  the  equivalent  circuits  (Fig.  4.12) 
the  following  equations  allow  us  a  unique  determination  of  the  required 
length  of  the  inductor  microstrip.  It  was  assumed  that 


2L}  +  L3  =  1.51  nH 
C  +  2CL  =  0.3  pF 
£  Z 

T  c  oc 

L1 

2  f  A 


f  A  Z 
gc  oc 


CL  2f  Z  T  A  T 
ol  gl 


A  ,  _  u)L~ 

£I  =  2n~  sin  ^ 

ol. 


Zqc  was  used  in  the  above  equations  (6)  as  an  unknown  parameter,  rather 

than  the  characteristic  impedance  of  the  5  mm  microstrip  betwewen  B  &  A. 
(At  A  &  B  there  are  microstrips  of  different  charat ter i sr  c  impedances  ) 
Next  the  width  'w'  of  the  inductor  microstrip  was  determined  from 


-  119.9 

/2j^r"TT) 


[£n{4  ^  +  V16(^)  +  2] 


- 1  <r+T><ta"  J  1  *4)] 


h  the  substrate  thickness  is  0.25  mm 
s 


Final  length  lj  and  the  width  w  of  the  inductor  microstrip  were  calcu¬ 
lated  to  be  1.1  mm  and  131  pm.  The  matched  transducer  was  then  fab. 
cated  through  usual  photolithographic  technique.  The  output  transducer 
was  a  mirror  impage  of  the  input  transducer  shown  in  Fig.  A.12After 
matching  the  following  improvements  have  been  observed 
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1)  It  has  been  possible  to  reduce  the  insertion  loss  by  6  dB  more, 

i.e.,  the  insertion  loss  at  resonance  with  most  of  the  1  mm 
resonators  can  be  kept  at  8-14  dB  over  3-6  GHz  (Fig.  4.16  At  3  GHz, 
insertion  loss  as  low  as  6  dB  at  resonance  has  been  obtained  for  some 
resonators  (Fig.  4.14) 

2)  The  out  off  resonance  rejection  has  been  improved  by  4  dB  i.e., 

10-14  dB  within  3-6  GHz  (Fig.  4.7,4.16,4.17).  Before  matching  it  was 
around  7-10  dB . 

3)  The  matched  parameters  have  been  found  to  be  quite  low  (0.15) 

around  5-6  GHz  than  compared  to  those  of  the  computer  optimized 
values.  In  the  3-5  GHz  range  measured  values  were  closed  to  the 
optimized  values. 

4)  Smith  Chart  plot  C  gives  the  measured  values  of  Sjj  with  one  of  the 

resonators  (identification  It  K2)  using  both  matched  and  unmatched 
1  mm  transducer. 

5)  Tunability  also  improved  after  matching. 

With  most  of  the  resonators  matched  transducer  did  not  provide  higher  Q 

resonance . 

Cone  1  us  ions  : 

1.  By  reducing  the  width  of  the  resonator  and  by  reducing  the  length  of 
loop  strip  of  the  transducer,  the  insertion  loss  at  resonance  is 
reduced  significantly.  This  is  also  reflected  in  the  lower 

parameters  (see  the  Smith  chart  plot  of  S.^). 

2.  By  matching  the  transducer  the  insertion  loss  has  been  reduced 
further  by  6  dB.  The  off  resonance  rejection  has  improved  by  at 
least  4  dB;  Q  however  did  not  improve. 

3.  Q  of  resonators  (as  studied  in  this  work)  depends  mainly  on  (i)  exact 
reproduction  of  the  mask  on  Y1G  i.e.,  uniform  land  and  Groove  spac- 
ings  throughout  the  arrays.  (ii)  large  II  of  bars  typically  ~  60 
(iii)  d/h  ~  2.5 %.  Earlier  workers  (Collins  et  aj ,  CASTERA  used  high 
d/h  (as  high  as  6-4%  to  11%)  to  obtain  high  Q.  Our  works  show  that 
high  d/h  basically  helps  to  obtain  a  single  moded  resonance  with 
large  off  resonance  rejection.  For  a  given  resonator  transducer 
combination,  Q  at  the  resonance  also  depends  on  the  position  of  the 
transducers  relative  to  the  spacings  between  the  two  arrays.  ey 
reducing  the  separation  between  the  two  arrays  the  spurious  response 
or  sidelobes  can  he  diminished  to  a  great  extent;  however  the  loaded 
Q  is  also  theoretically  predicted  to  be  reduced.  Reduction  in  Q  also 
occurs  with  the  reduction  in  the  width  of  the  resonator.  Reduction 
in  Q  may  be  compensated  by  properly  choosing  the  etch. 


Figure  4.14This  figure  shows  6  dB  insertion  loss  at  resonance. 

f=3.5  GHz;  1.5  mm  wide  resonator  (?iK  21)used  with  1  mm 
matched  transducer.  The  high  frequency  side  resonance 

however  was  not  observed  with  other  1  mi"  resonators.  9 


Figure  4.15  1  mm  wide  NSSW  resonator  (//  K2)  response  at  3  GHz  with  1  mm 
matched  loop  transducer.  The  insertion  loss  is  10  dB  at 

resonance.  The  bottom  trace  corresponds  to  ^ .  It  is  9 

seen  that  there  is  not  a  significant  dip  at  resonance.  . 

The  S  is  about  0.4.  The  Q  at  resonance  is  about  GOO. 

For  optimum  Q  we  have  to  .sacrifice  11.  by  2  dB.  The  '0 

corresponding  Q  at  G  GHz  was  1800  ar.d  11.  -  lb  dB. 

» 


depth  (d/h  ~  2-3%),  land  and  groove  width  as  well  as  the  position  of  the 

transducer  with  respect  to  the  cavity  spacings.  Our  works  show  that  low 

array  separation  (<  2  mm  does  not  provide  high  off  resonance  rejection). 

Higher  off  resonance  rejection  (~20  dB)  may  be  possible  by  carefully 

choosing  the  separation  of  the  loop  strips.  Because  of  metal  loading  by 

the  YIG,  the  effective  wavelength  A  of  the  wave  ,  ,  ,  , 

°  eff  launched  by  the 

transducer  differs  from  that  expected  on  the  basis  of  equation  Ao=2s 

where  s  is  the  separation  between  the  two  microstrips  of  a  loop.  If  this 

separation  can  be  chosen  carefully,  the  insertion  loss  will  improve 

further  and  will  be  less  than  14  dB  at  6  GHz  (Fig.  4.7)  This  is  expected 

to  improve  the  rejection  out  off  resonance,  too. 


It  may  be  mentioned  that  at  higher  bias  field  (which  corresponds  to 
higher  resonance  frequency,  the  transducer  passband  shrinks  because  of 
the  narrowing  down  the  difference  between  the  upper  and  lower  cut  off 
frequencies  in  ferrite  (see  Fig.  4. 3, 4. 4, 4. 5)  for  the  MSSW  wave. 

Because  of  this,  the  matching  of  wavelength  of  the  wave  launched  by  the 
transducer  with  the  periodicity  of  the  resonator  grating  becomes  very 
important.  A  small  mismatch  will  place  the  resonance  peak  at  a  different 
position  other  than  the  top  of  the  passband  and  this  will  add  to  the 
insertion  loss  significantly  at  high  frequency.  Change  of  wavelength 
because  of  the  metal  loading  by  the  YIG  needs  to  be  taken  into  account. 
These  above  mentioned  problems  are  found  to  be  very  significant  for  short 
narrow  resonators  (8  mm  long  and  1  mm  wide)  operating  at  low  wavelength 
(100  micron).  (Fig.  4.20) 


Experimentally  further  improvements  may  be  attempted  by  giving  ±  small 
variations  about  the  length  and  width  of  the  inductor  microstrip  (Fig.  4.12) 
to  reduce  the  at  high  frequencies. 

In  this  work,  single  modedness  of  the  resonance  over  3-6  GHz  was 
determined  by  finding  out  the  phase  shift  at  the  resonance  as  the 
resonator  is  tuned  over  3-6  GHz.  Only  phase  shift  about  10-12°  could  be 
recorded . 


Figure  4.16  1  nun  resonator  (MSSW)  response  (’inability  in  the  range 

2.75-6  GHz  with  a  1  mm  matched  loop  transducer  and  resonator 
#K2  is  used. 


Figure  4c 17  Resonance  response  ((inability  in  the  range  2.5  6  GHZ 
with  a  different  resonator  transducer  combination, 
wide  resonator  used  with  1  mm  matched  transducer. 
Resonator  <1 1  is  used. 


HSf 

u,  > 

Figure  4.18  l  mm  wide  resonator  response  with  1  mm  wide  matched 


transducer.  The  upper  curve  corresponds  to  trace. 


The  bottom  trace  corresponds  to  S 


The  center 


frequency  is  3.75  Ghz .  Sweep  range  =  1  GHz.  It  may  be 
noted  that  there  is  a  very  significant  dip  in 
occurring  at  the  resonance  point.  This  happens  only 
when  the  resonator  is  matched.  For  unmatched  resonator 
see  figure  below. 


Figure  4.19  1  mm  wide  resonator  (7/K2),  same  as  used  for  the  upper 

figure)  with  1  mm  wide  unmatched  transducer.  It  may  be 
noted  that,  the  significant  dip  in  Sjj  occurs  at 

different  frequency  than  the  resonance  point.  Resonator 


rating  at  4.o  ohz. 

,  i gned  for  X  =  220  pm. 

>a ,  the  resonance  is  pi  =» «' 
1.  The  transducers 
g  made  on  the  off 


]  mm  x  o  uun  x 
At -1  GHz .  The 
It  is  seen  tha 
at  the  top  of 
are  unmatched, 
resoname'  rt\)f. 


Single  Mode  Character  of  Resonance 

In  this  work  single  mode  character  of  the  resonance  was  determined  by 
finding  out  the  phase  shift  relative  to  the  resonance  at  3  GHz  as  the 
resonance  frequency  is  varied  to  6  GHz  at  steps  of  200  MHz.  When 
corrected  for  the  transmission  phase  (^2^)  (of  the  signal  passing 

through  the  two  matched  transducer  feeder  lines  (Fig.  4. 12  as  obtained 
from  compact  (see  Fig.  lib')  we  find  that  there  is  only  overall  phase 
shift  of  10-20°  as  the  resonance  is  tuned  over  3-6  GHz  by  adjusting  the 
magnetic  field.  (See  Fig.  4.22) 


This  small  phase  shift  might  be  due  to  the  fact  the  compact  did  not  take 
into  account  of  the  frequency  dependent  electrical  characteristics  of  the 
microstrips.  Further  considering  the  error  in  the  measurement  of  phase 
(which  is  typically  around  10°)  we  can  say  that  the  resonance  is  single 
moded . 


Fiq  11a Phase  of  S„,  measurements  with  resonator# K2(l  inn  wide)  at  resonance. 

1  m  matchii  transducer  Is  used.Callbratlon  Is  done  with  microwave  short. 

xxx  phase  at  resonance  Including  that  of  transducer  feeder  lines, 
ooo  phase(electr1cal  lenght  of  two  Identical  cascaded  feeder  lines. 
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Figure  4.22  Phase  of  S,,  measures  with  resonator  #K2  (1  mm  wide) 

at  resonance.  1  nun  matched  transducer  is  used. 
Calibration  is  done  with  microwave  short. 

xxx  phase  at  resonance  including  that  of  transducer 
feeder  lines. 

ooo  phase  (electrical  length  of  two  identical  cascaded 

(nA/lar  1  i  nPC 


Electrical  length  of  the  input  and  output  cascaded 
Feeder  lines  of  the  transducers. 

BEDIT  (FILE)  KK2.DAT.3  (OUTPUT  PS)  KK2.DAT 
Edit)  KK2.DAT.3 

*P 

80100  MET  AA  ZZ 

00200  TRL  AA  MS  ,52  5  3.8  .254 

00300  CAP  BB  PA  -.30114 

00400  IND  CC  SE  -1.5096 

00500  TRL  DD  MS  .55  6  3.8  .254 

00600  CAX  AA  DD 

00700  CAS  AA  AA 

00800  PRI  AA  SI  50 

00900  END 

01000  3000  4400  200 

01100  4600  4000  5000  5100  5200  5300  5400  5500 
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01300  END 
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Figure  4.23 


Magnitude  and  phase  angle  of 

feeder  lines  of  the  1  matched 
Page  170) 


^21  ^or  ^wo  cascaded 
transducer  (see  fig.  4.12 
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COMPACT  matched  S,.  parameters  for  Imt  resonator.I  studied  with 
1  rat  X  ducer  1 
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I  NAME 


'TITLE 


Smith  Chart  plot  of  measured  at  resonance  point  over  3-6  GHz.  1  mm 
resonator  #1  was  used  with  1  mm  unmatched  loop  transducer.  Matching 
through  COMPACT  was  done  on  these  parameters.  Frequency  increases 
alphabetically.  (See  page  190)  Corresponding  transducer  short  was  used 
for  calibration. 
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Resonator  Delay 

It  may  be  mentioned  that  we  have  carried  out  time  delay  measurements  at 
resonance  by  finding  out  the  phase  shift  when  the  frequency  is  changed 
by  Af  (=  5  MHz)  about  the  resonance  point.  The  delay  is  measured  through 
the  equation 

r  =  _L_ 

1  2ttL  Af 

where  L  is  the  spacing  between  the  input  and  output  transducers.  The 
measured  time  delay  over  the  frequency  range  3-6  GHz  is  shown  in  Fig.  for 
the  resonator  #1  and  a  matched  transducer  (1  mm  microstrip).  THe 
resonator  transducer  combination  was  set  for  Q  =  550  -  1700  in  the  range 
3-6  GHz  with  corresponding  insertion  loss  8-16  dB.  Increase  of  insertion 
loss  by  2  dB  (higher  than  14  dB  obtained  for  fresh  transducer)  is 
possibly  due  to  the  deterioration  of  the  transducer  microstrip  and  the 
OSM  connectors. 


It  is  found  that  the  delay  varies  linearly  (Fig.  4.24with  resonance 
frequency  in  the  range  3-6  GHz.  It  should  be  noted  that  the  different 
delay  times  corresponds  to  different  bias  fields. 
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Fiq.4.24  T1nie  delay  measurements  at  resonance 

Resonator# I ( 1  mn  wide)  and  1  mn  matched 
transducer  used. 
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Fig.  4.24  Time  delay  measurements  at  resonance;  Resonator  #1  (  1mm 
wide)  and  1  mm  matched  transducer  used. 


4.1.4 


Physical  Characteristics  of  the  Transducers 


1.  Resonators  with  3  mm  aperture  are  studied  with  3  mm  long  x  ducer 

(i.e.,  the  loop  microstrips  are  3  mm  long).  The  separation  between 
the  input  and  the  output  loop  strips  is  2.4  mm.  The  microstrip 
feeder  lines  are  530  micron  wide.  The  total  resistance  between 
input  and  output  OSM  connectors  was  4.8  ohm. 

2.  Resonators  with  1.5  mm  aperture  are  studied  with  1.5  mm  long  x 

ducer.  The  separation  between  input  and  output  loops  strips  is  2.1 
mm.  The  total  resistance  between  input  and  output  OSM  connectors 
was  1.6  ohm . 

3.  The  1  mm  wide  resonators  were  studied  with  1  mm  long  x  ducers.  The 

separation  between  the  loop  strips  is  2.4  mm.  The  total  resistance 
was  0.9  ohm. 

All  the  x  ducers  are  fabricated  on  quartz  substrates  of  thickness 
0.254  mm  and  made  of  gold  microstrips. 
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Physical  Characteristics  of  the  Resonators 


Summary  of  Resonator  Response  Data 


Summary  of  Resonator  Response  Data  (Cont) 


4.1.5 


ABSTRACT 


YIG  Oscillators:  Is  Planar  Geometry  Better 

by 

R.  L.  Carter 
J.  M.  Owens 
D.  K.  De 

With  the  advent  of  planar  magnetostatic  wave  (MSW)  resonators  there  are 
two  YIG  resonator  technologies.  The  planar  technologies  will  be  compared 
to  the  YIG  sphere  technology  with  regard  to  design  and  fabrication 
consideration.  New  Data  on  the  optimization  of  (MSW)  resonators  will  be 
given. 


Broad-band  tunable  signal  sources  at  microwave  frequencies  are  an 
important  systems  requirement.  In  particular,  communications  systems 
require  oscillators  with  high  spectral  purity  and  low  phase  noise.  A  new 
planar  yttrium  iron  garnet  (YIG)  technology  has  been  developed  over  the 
last  ten  years  which  utilizes  the  propagation  of  magnetostatic  waves 
(MSW)  in  thin  expitaxial  films  of  YIG  on  gadolinium  gallium  garnet 
(2) 

substrates.  Recently,  we  have  built  resonator  structures  in  the 
YIG/MSW  technoigoy  which  have  been  used  in  oscillators  reported  to  have 
phase  noise  characteristics  of  -100  dB  at  10  kHz  offset  from  the 

Ill  C 

carrier  frequency.  Castera,  et  al ■  have  reported  -95  dBc .  Tunable 

MSW  oscillators  have  been  reported  form  2  GHZ  to  10  GHz  with  multiple 
octave  tunability.  For  more  than  30  years,  YIG  based  resonators  have 
been  used  in  signal  sources  utilizing  spherical  YIG  structures.  This 
paper  will  compare  the  strengths  and  weaknesses  of  the  application  of 
planar  versus  spherical  structures  in  the  design  and  fabrication  of 
micro-wave  oscillators. 

The  consideration  of  factors  of  material  related  parameter  il 
apply  similarly  lo  both  the  planar  and  sphere  geometry.  Howevi  with 
regard  to  phase  noise  considerations,  the  variance  is  apparent  sine"  the 
structure  determines  the  resonator  characteristics.  The  YIG  sphi  •' 
resonator  is  responding  at  the  uniform  spin  precession  frequency.  The 
planar  YIG/MSW  resonator  is  propagating  energy  at  a  MSW  mode  frequency. 

As  a  consequence,  the  YIG/MSW  resonator  will  always  have  a  higher 
resonant  frequency  for  a  given  magnetic  field  than  the  YIG  sphere 
resonator.  The  most  crucial  difference  appears  in  the  phase  noise 
characteristics.  For  a  resonator  stabilized  oscillator,  the  phase  noise 

(3) 

can  be  modelled  by. 

G  k  T  F  F 

L(f )  =  10  log  -x - -^5 

Q  P  (f-f  ) 
o  o 


where  G  is  the  amplifier  gain  (or  IL,  the  insertion  loss  of  the 
resonator),  kT  is  -174  dB^,  F  is  the  noise  figure  of  the  amplifier,  fQ  is 
the  resonant  frequency,  f  is  the  frequency  at  which  the  noise  power  is 
measured,  Q  is  the  resonator  quality  factor  and  Pq  is  the  output  power 

level  of  the  amplifier.  In  comparing  two  saturated  resonator  based 

2 

oscillator  technologies,  the  factor  IL/Qg  Pg  is  the  figure  of  merit, 
where  Qg  is  now  the  saturated  resonator  quality  factor  and  Pg  is  the 
resonator  saturation  power.  Of  particular  interest  is  the  factor  Qg : 

The  reported  saturated  quality  factor  of  the  planar  YIG/MSW  resonators 
exceeds  that  of  YIG  spheres  by  nearly  a  factor  of  10.  The  paper 
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Figure  4.25  Oscillator  block  diagram 
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will  discuss  the  limits  of  each  technology  with  regard  to  IL,  Q  and  P  , 
with  new  data  being  reported  on  methods  of  the  optimization  of  ?L  and  § 
for  YIG/MSW  planar  resonator. 

The  fabrication  techniques  required  for  an  oscillator  technology  are 
of  considerable  importance.  The  paper  will  discuss  the  following  factors 
relating  to  the  planar  versus  spherical  technology. 

1.  Process  Yield 

2.  Interface  to  the  Oscillator  -  is  integration  possible? 

3.  "Tweaking”  processes 

Some  of  the  MSSW  resonators  fabricated  in  this  work  were  driven  by  a 
GaAs  FET  amplifier  to  oscillate  between  3-5.3  GHz.  The  hybrid  integrated 
amplifier  was  fabricated  on  10  m£  alumina  substrates  with  gold 
metalization.  The  resonator  and  hybrid  MIC  amplifier  were  assembled  into 
the  oscillator  configuration  of  Fig.  4.?Fi  The  two  stage  amplifier  had  a 
gain  of  18-19  dB  in  the  range  3-5  GHZ. 


Phase  noise  characteristics  of  some  of  the  oscillations  are  shown  in  Fig. 
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figure  4.26  The  above  oscillations  wore  meaasured  with  a  combination 
of  RESN  4i K2  and  1  mm  matched  Xdncer  built  on  quartz 
substrate  of  thickness  250  micron.  The  above 
combination  was  placed  in  feed  back  loop  of  a  two  stage 
microwave  microstrip  amplifier  fabricated  in  this 
laboratory  on  250  micron  thick  alumina  substrate.  The 
amp]  i  tier  has  a  gain  of  18-19  dB  in  the  range  of  '1-5  GHz 


5.0  SUMMARY 


In  summary  the  primary  accomplishments  of  the  third  year  of  this 
program  have  been: 

1.  An  Empirical  Model  has  been  developed  for  MSSW  Transducers  using 
detailed  measurements  of  MSSW  transducers.  The  theory  is  still  limited  by 
errors  in  measurement  and  the  inability  to  measure  S^  and  S^2  of  the 
magnetostatic  ports.  Also,  errors  associated  with  phase  shifts  under  the 
transducers,  and  interaction  between  the  transducer  elements  are  not 
considered.  The  solution  to  the  problem  lies  in  either  detailed 
theoretical  analysis  of  the  problem  (6  region  boundary  valve  problem)  or 
decoupling  of  the  transducer  elements  to  yield  a  true  non-recursive 
transducer  filter. 

2.  A  variable  time  delay  device  with  an  RMS  phase  error  of  as  low  as 
6.40°  RMS  has  been  built  and  tested.  Time  variation  of  42  NSEC  was 
achieved.  The  design  procedures  for  dispersion  control  utilizing  a 
variable  ground  plane  are  well  established.  The  device  utilizes  a  cascaded 
forward  and  backward  wave  delay  line  with  linearized  dispersion,  and  thus 
matched  dispersive  line  have  been  developed  with  significant  potential  in 
signal  processing  applications. 

3.  Two-port  MSSW  resonators  have  been  developed  with  loaded  Qs  of 
greater  than  2500  at  6  GHz.  Oscillators  built  from  these  resonators  has 
shown  excellent  phase  noise  characteristics  with  SSB  phase  noise  as  low  as 
-  105dBC  at  10  KHz  offset  from  carrier  and  -135dC  at  1  MHz  offset. 
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